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Abstract 

We apply our previous work on Green's functions for the four-dimensional quaternionic Taub-NUT 
manifold to obtain a scalar two-point function on the homogeneously squashed three-sphere (otherwise 
known as the Berger sphere), which lies at its conformal infinity. Using basic notions from conformal 
geometry and the theory of boundary value problems, in particular the Dirichlet-to-Robin operator, we 
establish that our two-point correlation function is conformally invariant and corresponds to a boundary 
operator of conformal dimension one. It is plausible that the methods we use could have more general 
applications in an AdS/CFT context. 
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1 Introduction 

The AdS/CFT correspondence |21 [3] provides a very precise example of the holographic phenomenon, 
whereby a d + 1-dimensional gravitational theory has an equivalent description as a d-dimensional theory 
without gravity. In AdS/CFT this "dual" theory is actually a gauge theory. In the limit where on one 
side of the duality we have weakly-coupled gravity the gauge theory on the other side is expected to be 
strongly-coupled and to have a large number N of colours. Essentially, as first observed in 0, in this limit 
one is simply dealing with a classical boundary problem for fields (perhaps with nonzero mass or spin) in 
AdS, whose (appropriately rescaled) boundary values are interpreted as sources for the dual fields in the 
gauge theory. The gauge theory can thus be thought of as living on the boundary of AdS. The bulk fields 
obey their equations of motion, so they can be constructed using Green's functions of the appropriate bulk 
laplacian. This allows the calculation of correlation functions in the boundary theory. We refer to 0] EI for 
reviews and a survey of the literature. 

The Einstein metric on Euclidean Anti-de-Sitter space, or hyperbolic space, is the prototypical example 
of a conformally compact Einstein metric on a compact riemannian manifold with boundary, defined as a 
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complete metric in the interior of the manifold which blows up with a double pole at the boundary, and 
can thus be extended to a boundary metric via a suitable choice of defining function^ . A defining function 
is any function p which is positive in the bulk, but has a simple zero at the boundary, so that p^g extends 
smoothly to the boundary. Since any such function can play the role of a defining function, the induced 
metric at the boundary is not unique. Thus we obtain an induced boundary conformal structure rather 
than a boundary metric. This is of course crucial for AdS/CFT to work, since any theory defined on the 
boundary will automatically be conformally invariant (modulo possible conformal anomalies in even boundary 
dimensions). However, this fact is true for any theory conjectured to be dual to gravity on a conformally 
compact manifold, so given the relative scarcity of conformal field theories in more than two dimensions, it 
is of interest to see how AdS/CFT can be extended to this setting and possibly provide nontrivial examples 
of CFT's. A first step in this direction would clearly involve solving the bulk laplacian to obtain Green's 
functions, which would then lead to CFT correlation functions, in a similar way to the AdS case. 

One way to approach this problem is the method of holographic renormalisation, reviewed in [7], which 
uses extensively the framework of Fefferman and Graham in which the bulk metric is expanded in a 
series of the radial distance from the boundary. One can do the same for any bulk field, by prescribing a 
source field on the boundary and then performing an iterative near-boundary analysis, in which the bulk 
field equations are solved algebraically and substituted back, thus providing higher terms in the radial series. 
Holographic renormalisation is a very general method, applicable to any asymptotically hyperbolic manifold, 
and is ideally suited to calculating the local counterterms needed to regularise and renormalise the divergent 
bulk action. However, its usefulness in calculating correlation functions in the boundary is limited, because 
the necessary coefficient in the radial expansion, being nonlocally related to the prescribed source, cannot 
be fully determined by the near-boundary analysis. At that stage the analysis needs to be supplemented by 
an exact solution of the bulk field equations for the field in question. So although one could go rather far in 
analysing the dual field theory using holographic renormalisation, it is still important to obtain exact bulk 
Green's functions for asymptotically hyperbolic manifolds. 

In in], we discussed a particular example of such an asymptotically hyperbolic manifold, the Quater- 
nionic Taub-NUT space (QTN). It is a four-dimensional manifold with SU(2) x U(l) isometry group, whose 
boundary conformal structure corresponds to a Berger sphere, which is a homogeneous squashing of the 
three-sphere preserving an SU(2) x U(l) isometry group out of the round three-sphere's S0(4). QTN can 
be thought of as a deformation of hyperbolic space HH^ (Euclidean AdS4), but a very special deformation 
in that it is conformally self-dual, meaning that the Weyl tensor (zero for HH^) is self-dual. Although a 
non-zero Weyl tensor means that, unlike the case of HH^ and its boundary the round sphere, the boundary 
manifold is not guaranteed to be conformally flat (and indeed it isn't), and thus deprives us of the standard 
trick of supposing the boundary is really just flat space, keeping the Weyl tensor self-dual is extremely useful 
in that there is a twistor description of the manifold. This means that we can translate differential-geometric 
questions about the four-manifold (in this case, finding Green's functions) into complex-analytic questions 
about its 6-dimensional twistor space, hoping that the more powerful techniques thus available will allow us 
to answer them. 

The outcome of ^ was the explicit calculation of bulk-to-bulk and bulk-to-boundary propagators for a 
conformally coupled scalar on QTN. Assuming there is a holographically dual field theory, we can apply the 
holographic prescription |2| to immediately find two-point correlation functions of the dual gauge-invariant 
operator, which (as will be discussed) is expected to have (mass) conformal dimension A = 1 in the strongly- 
coupled boundary theory. We will write down the two-point function in section |21 Checking that this 
two-point function, which we call M{x,x'), is indeed the "correct" one (i.e. verifying that it is conformally 
invariant, and corresponds to a A = 1 operator) will occupy the rest of this article. 

Why is it not immediately obvious what conformal dimension our correlator corresponds to? After all, 
in flat space one only has to look at the power of the l/jx — x'| term in a scalar two-point function and read 
off the conformal dimension of the corresponding operator. However, recall that all this information arises 
through the analysis of representations of the conformal group. As is well known, on the round sphere (or, 
equivalently, fiat space) this group is large enough to allow a complete characterisation of two-point and 
three-point correlation functions (see e.g. ^Hl)- As an illustration, suppose M{x,x') is a two-point function 

^The condition of conformal compactness, together with the Einstein condition, imphes that the metric is asymptotically 
hyperbolic, that is, all sectional curvatures approach —1 near the boundary 10. (For hyperbolic space, the sectional curvature 
is —1 everywhere). So metrics of this type are often also called asymptotically hyperbolic, or asymptotically AdS. 
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(in d dimensions) of the scalar operator O, which has conformal dimension A. Then, if u is a conformal 
Killing vector, we have the following condition on the correlation function^: 

+A(a„(a;) +a„(a;'))] M{x,x') = (1.1) 

where Cv is the Lie derivative, and a^ix) is defined through the conformal Killing equation 

jO-v{x)9fj.u - '2<7vix)gfj,,y ^ (1.2) 

to be cr^ = V • f/d. As we will discuss in section O if we denote by s the order of the differential operator 
T>x associated to M{x,x'), in the sense that 

V,M{x,x') = ^S^^\x,x') (1.3) 

(in other words, V is the inverse of M(x,x')) then the quantities d, A and s are related by 

2X^d-s (1.4) 

As an illustration of (|1.1() . pick a flat-space conformal Killing vector, for instance the dilatation vector 
V = a;^(9^, and apply it to the well-known flat-space conformal correlation function M{x, x') ^ l/[(a;— x')^]^. 
We get (notice here crv{x) = cr^(x') = 1) 

[t^v{x) + l^v{x')\M[x, x') = -2AAf(x, x) (1.5) 

corresponding to the well-known fact that M(x, x') has conformal dimension A. If we put s = 2, i.e. M{x, x') 
is the Green's function of the laplacian, we get the dimension A — d/2 — 1 of a free scalar field. 

Now let us turn to the Berger sphere. By Obata's theorem VZ^ the only compact manifolds where 
the conformal group is larger than the isometry group are the spheres, equipped with the standard bi- 
invariant metric. So in the case of the Berger sphere, which is only left-invariant, the conformal group is 
just the isometry group, in other words there are no conformal Killing vectors apart from the Killing vectors 
themselves"^. So the method just outlined would not work for the Berger sphere, and we need to look for 
more sophisticated ways of showing that our correlation function has the form that is required by conformal 
invariance on the Berger sphere. 

It turns out that a suitable framework to discuss this problem is that of conformal geometry, a field 
pursued in recent years by several mathematicians, and which has deep connections to the ideas of Fefferman 
and Graham on conformal invariants. Thus section|31will be devoted to a lightning review of the main ideas. 
Then, motivated by H1.3|l . we will be led to ask the question: "Since our correlator is a Green's function, which 
is the operator it is the Green's function of?". In other words, given our two-point function M{x,x'), we 
will look for an operator T) on the Berger sphere such that H1.3(l holds. The idea is that, although Hl.l|l is not 
useful anymore in constraining M{x, x'), if we can find an operator V, with known conformal transformation 
properties, such that equation (|1.3|l holds, then it can be used to constrain M(x,x'). Of course this hinges 
on the existence of an independent way of finding V. Note that if M{x,x') indeed corresponds to A = 1, 
then from H1.4(l we should expect I? to be a first-order operator. 

In section 0] we will see that this "D turns out to be the Dirichlet-to-conformal-Neumann operator (we 
will mostly call it Dirichlet-to-Robin for short), which for hyperbolic space takes the simple form: 

B - ^Ab + 1/4 (1.6) 

where is the (positive) laplacian on the boundary. So S is a pseudodifferential operator of order one on 
round S^, as anticipated'*. Before explaining why, on general principles, B is expected to be the inverse of 

^We will very briefly indicate how this equation comes about in section |21 A detailed discussion of similar formulas can be 
found in 1111 . 

^Note however that in the limit of extreme prolate squashing the left-invariant sphere degenerates to a CR structure, which is 
not a smooth manifold, Obata's theorem does not apply, and indeed (as discussed in 1131 . where it was realised as the boundary 
of the Bergman space CH'^) this space allows conformal Killing vectors. Lee 1141 and Schoen 1151 have analysed this case (the 
unit sphere in C"), which turns out to be the unique CR manifold to admit a nontrivial conformal group. The analogue of 
conformal flatness in this case is the CR equivalence to the Heisenberg group. 

*It should be clear that B has nothing to do with the Dirac operator. Apart from the fact that = A5 + I/4 = 
Ab + TZb/Q (the scalar curvature of S"^ is TZb = 3/2 in our conventions — see appendix while (according to Lichnerowicz) 
= Ab + 7?._b/4, it is crucial that B has to be a nonlocal operator, as will be discussed in section |3 
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M(x,x'), we will discuss a way to construct this operator for a given manifold with boundary, give some 
background on the kind of problems it appears in, and go over its properties as a conformally covariant 
operator. 

The outcome of section0]is that we can verify that M{x, x') is conformally invariant by simply checking 
with T) = B. We find it instructive to demonstrate this first in the simple case of hyperbolic space 
where everything is clear and can be easily understood. This is the purpose of section [S] To act with B on 
M{x,x') we will find it convenient to expand M{x,x') in eigenfunctions, which is trivial to do on the round 
sphere. 

The confirmation we have been looking for will finally appear in sectional where we apply our methods 
to Quaternionic Taub-NUT, and its conformal boundary the Berger sphere. Here too the fact that the 
eigenfunctions of the laplacian are known will be crucial in avoiding the use of heavy functional-analytic 
machinery. The drawback is that expanding M(x,x') is a rather more involved technical exercise. 

Section 13 contains some conclusions and possibilities for future work. We have also included an appendix 
collecting various facts on the squashed three-sphere which we will use at various places in this paper, as 
well as an appendix with some facts and references on pseudodifferential operators. 

2 The conformally invariant correlation function 

In this section, after first discussing some properties of the QTN metric, we state the results of the analysis 
in [HI, namely the bulk-to-bulk and bulk-to-boundary propagators for a conformally coupled scalar. Then 
we can easily find a two-point function on the boundary Berger sphere, which we argue should correspond 
to a dual operator of conformal dimension A = 1. We take the AdS limit and check that it is so in that 
special case. 

2.1 The Pedersen metric for QTN 

Our starting point is the Pedersen metric for quaternionic Taub-NUT ^: 

_^2fc2^4 d^ +r [l-^ir )ia,+a^) + ^^^2^2 '^aj ■ (2-1) 

The parameter jj? can take values in the range [k'^, —00), while the parameter fc^ is related to the (negative) 
cosmological constant through A — — 3fc^. Keeping jj? e [0, fc^] corresponds to a prolate squashing of the 
three-sphere at the boundary, while jj? E [0, —00) to an oblate squashing. We will restrict our attention 
to the prolate case in this article. In this range QTN has two limits which are special, in the sense that 
they are symmetric spaces: /i^ = is just 4-dimensional hyperbolic space (in other words, Euclidean AdS4) 
which we denote by HH^, and /i^ = is an Einstein-Kahlcr manifold known as the Bergman space, often 
denoted by CH^, a notation we will adopt^. As cosets, these manifolds are HH^ = S0(1, 4)/SO(4) and 
CH2 = SU(2,1)/U(2). 

As for the radial coordinate r, we take it to lie in the range [0, 1/fc). This choice ensures that the metric 
is complete. The boundary is at r — > 1/fc where the metric blows up and smoothly induces a boundary 
conformal structure. Finally, the ai are the left-invariant SU(2) one-forms, which are given explicitly in 
appendix IXI Taking (1 — fc^r^) as our defining function, and calling i,j the non-radial directions, the induced 
metric at the boundary is 

= ^lim^(l - fcVf g., d.2 = l[(l-MVfc^)(a? + a|) + cT|] (2.2) 

This is the left-invariant metric of the Berger sphere, which is examined in some detail in appendix 

The twistor space of QTN was constructed in 16 . Knowledge of the twistor space was essential in P], 
since it allowed us to make an educated ansatz (following a similar construction of Page ,17 for Euclidean 

^There are several other forms of the metric for this manifold, which is often called AdS-Taub— NUT in the physics literature. 
See |5| for references and more details. 

Sjn the Ber gman space was denoted by CP^. 



(l-fc2r2)2 
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Taub-NUT) for the Green's function in the case of a conformally coupled scalar field. Using this ansatz we 
were led to an explicit form for the Green's function. In this article we will start with this result, so no 
twistor methods will be used. 

We now wish to consider free scalar fields propagating on the QTN geometry^, having a conformal 
coupling to the curvature. We recall from jSj the scalar laplacian corresponding to the metric H2.1|l : 

V9 



(1 - n^r^) 



(1 - M'fcV^) dr 



4r(l - k^r^) 



(3 + k^r^ - 7^?k^r^ + S^^/cV^) dr 



(2.3) 



cot Ode + csc^ 9 



d^^ - 2 cos 68^^, 



sm 



1 — fi^k'^r'^ 



cos 



The conformal laplacian Y ~ — TZ/6 corresponds, since the scalar curvature is 7?. = 4A = — 12fc^, to the 
operator Y = V'^ + 2k^. We should now supplement this operator with appropriate boundary conditions on 
the boundary of QTN. 

Note that for massless scalars, with an arbitrary curvature coupling ^, satisfying Dirichlet boundary 
conditions, the presence of the bulk TZ term leads us to consider the following action (see [T8]): 



S 



JM JdM 



(2.4) 



where /C = — VriV'" is the trace of the extrinsic curvature (iV is the outward unit normal vector) and Ziy- 
is the (unrescaled, thus divergent) restriction of the metric to the boundary. The term with the extrinsic 
curvature plays a similar role to the Gibbons-Hawking term |19l in the Einstein-Hilbert action, i.e. it ensures 
that the variational problem is well defined, because with this addition: 



1. 5 is stationary under g^^^ 
boundary, and 

2. 5' is also stationary under c, 



diiv + 5g^^, with the metric and its normal derivatives fixed on the 

■ (j) + <^0j where (j) satisfies Dirichlet boundary conditions 5(j)\QM ~ 0- 

Now recall that, as first discussed in (see also ^1|^), for low values of the conformal dimension we 
have a choice of boundary conditions in AdS/CFT, in the sense that the same solution for the bulk field can 
correspond to two different boundary operators, depending on the asymptotic behaviour we wish to impose. 
The "regular" boundary condition (which can be imposed for any bulk mass) leads to a boundary operator 
with the usual conformal dimension A-|_ = d/2 + ^/dF/A + rn^, while the "irregular" boundary condition 
leads to A_ = d/2 — \J <P jA + vn?-. Massless, conformally coupled fields (where effectively — SJZ, and TZ 
is negative for asymptotically hyperbolic spaces) lie in this region of two possible quantisations, however one 
should check whether both boundary conditions are consistent with the conformal coupling. Since unitarity 
in a CFT requires A > (d — l)/2, and also the fields in the bulk should have masses above the Breitenlohner- 
Freedman bound of — d^/4, it follows that the irregular modes can only exist in the range < d^ /A + SJZ < 1. 
However, the analysis of the irregular modes in ^Sj leads to a further restriction on the values of ^: ^ = 
where A_ = d/2 — ^ d^ /A + SJZ. For massless scalars, there are two solutions for ^ in d + 1 = 4 
(which is not interesting since A_ takes the non-unitary value of zero) and C = ^ (which gives A_ = 
the conformal value ^ = 1/6 is allowed, and thus it is consistent to impose Dirichlet boundary conditions for 
conformal coupling. 

On the other hand, as dicussed e.g. in (55!rT5^, it is a consequence of the conformal coupling that we 
cannot consistently impose Neumann boundary conditions for the action H2.4|l . Varying the action with 
respect to (j) (and using the bulk equations of motion) we have 



2d 

1). So 



LdM 



(2.5) 



'^That is, we ignore self— interactions as well as any backreaction on the geometry. 
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So, while we are obviously free to impose the Dirichlet condition, the only other boundary condition that 
keeps the action invariant is the Robin-type condition: 

(^iVa,. -i/C^ ^Oasr->l/fc. (2.6) 

Here we have already substituted the conformal value of ^. The precise scaling as we approach the boundary 
will be dicussed below. In what follows we will sometimes call this condition conformal Neumann, but will 
frequently revert to calling it simply "Robin" for short. 



2.2 Green's functions on QTN 

Before writing down the Green's functions on QTN constructed in [3], we recall a few basic facts about 
boundary problems and Green's functions in the asymptotically hyperbolic case. This is well known material 
in AdS/CFT, but since the setting here is slightly more general we find it useful to retrace some steps. 
Given an asymptotically hyperbolic manifold Ai in d + 1 dimensions, the massive laplacian — ~ 
+ X{d — X)k^ admits two independent solutions which, in coordinates such that the boundary is at y = 0, 
scale as and y'^~^ (see e.g. [23 )■ So the general solution for a massive scalar behaves asymptotically as 

0(y, x) = a?) + y^Hy, x) . (2.7) 

The first terms /o ~ /(O, x) and Hq — h{0, x) in an expansion of /(y, x) and h(jj, x) are the only independent 
quantities, with the higher orders being uniquely (and locally) specified by them 120117]. Now for large values 
of the conformal dimension (meaning A > d), the first term in (|2.7(l is non-normalisable and corresponds to 
(modified) Dirichlet boundary conditions, while the second one is normalisable and corresponds to (modified) 
Neumann boundary conditions. So the Dirichlet problem consists of specifying /o on the boundary, while 
in the Neumann problem we specify /iq. In each case we can construct the solution for (j){y,x) in terms of 
a suitable Green's function. As discussed in |22I the asymptotic behaviour of these Dirichlet and Neumann 
Green's functions is — > and G'^ — > y'*"'^ respectively. The asymptotic behaviour of the Neumann 
Green's function comes from the requirement that it satisfy {N^dy + (d — \)k)G\ — > y^ asymptotically. 

Now we turn to our case of a conformally coupled scalar in d + 1 =4, where effectively m? — — 2fc^, 
resulting in the two solutions A+ = 2 and A_ = 1. Let us then write the asymptotic behaviour of (/)(y, x) as 

c^{y,x)^yf[x)+y^h[x) + --- (2.8) 

We have dropped the subscript of / and ft,, which should cause no confusion since the higher order terms 
in each are contained in the "• • •". As discussed in the previous section, in this case, where both modes 
are normalisable, we have a choice of boundary value problems, depending on which A we pick. Suppose we 
want to solve for (/)(y, x) by specifying h{x). If A = 2, this is clearly a Neumann problem. However, thinking 
of it from a A = 1 point of view, it is the Dirichlet problem. We recall that the (modified) Dirichlet problem 
in this case is the question of constructing a function 0(y, x) satisfying 

V^ + 2e)^{y,x)^Q 
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while the (modified) Neumann problem corresponds to finding a (t>{y, x) such that' 

V2 + 2P)</,(y,f)=0 

y-\N + (d - A)fc)0(y, x)] = ~h{x) ^^'^^^ 

where N — N^dy is the outward unit normal, as before^ . By [• • ■]dM we mean the restriction of the point 
y to the boundary, and [• • ■]dM' below will mean the restriction of y' to the boundary. Finally, if we know 

^The extra term in the boundary condition can be easily seen to follow from Green's third identity applied to the modified 
Neumann problem 1221 . 

^It is basically given by —ydy in the {y,x) coordinates. 
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the Dirichlet and Neumann Green's functions, the Dirichlet problem can be solved through a double-layer 
potential (here dQ{x) is the measure of the induced, finite boundary metric) 



0(2/,^) 



[{y')-\Ny, + id- X)k)G^o (y, f , y', x')] • f{x')d^l{x') 



while the Neumann problem is solved using a single-layer potential 



[(y')^-''G^^{y,x,y\x')]aM' ■ h{x')dn{x') 



(2.11) 



(2.12) 



As we will see explicitly in sections l5. II and l5. II for A = 2 the factor of (d — X)k is precisely the restriction 
to the boundary of the extrinsic curvature term — 1/3/C appearing in the Robin problem, so we can actually 
switch to calling Gff{y,x,y' ,x') the Robin Green's function, and denote it Gj^{y,x,y' ,x'). 

After these preliminaries, let us turn to the bulk-to-bulk Green's function for a conformally coupled 
scalar on QTN, which was found in ^ to be^° 



Gi^^^Hr,x,s,x')^^'-''^'^^--''''^ 



Q+~Q- 



}+ + Q-- 4rsCvi 



where 



(2.13) 
(2.14) 

(2.15) 



with D = (r^ - s2)2 - 4r2s2(v2 - 1) + AiPk^r'^s'^Y2{y2 - 1) and 

G=[{1- ^i^k^r^){l - M'fc's*)] v/(l + Mr-2)(l + /xfcs2) . 

The vi and V2 are SU(2) x U(l)-invariant combinations of the three Eulcr angles describing the three-sphere 
(about which see appendix 

9 6 1 9 9 1 

Vl = cos y cos y cos -(^Ar + - '>Ps - Vs) + sin y sin y COS - (^^r " '/'r ^ V's + Vs) (2-16) 



and 



V2 



-{1 + cos9r cos9s + sm9r sm9s cos{ipr — ^s)} 



(2.17) 



Note that we have ignored the normalisation of H2.13I) and we will similarly ignore normalisation constants in 
the following propagators. G^^'''''> satisfies y/g{W^ +2k^)G'^'''''^{r,x,s,x') ^ S^'^^Xr^Xs) in the bulk of QTN, 
and it scales as (1 — k^r^) at the boundary r 1/k. So, according to the discussion above, if we have A = 1 
it is the Dirichlet Green's function Gjj, while for A = 2 we should think of it as the Robin Green's function 
G|j.. We can also verify for A = 2 that its Robin derivative scales as 



{N + k)G^''-''\r,x,s,x') 



(l-fcV)2 



dM 



(2.18) 



as it should. Now we can do two things starting from G^^'''^{r,x,s,x'): Thinking in A = 1 terms, we can 
solve the Dirichlet problem by writing 



r, .T = - 



dM' 



(l-fc2s2) 



{N''d, + 2k)G^'''''\r,x,s,x') 



Mx')dn{x') 



(2.19) 



dM' 



Note that 4>o{x) is really h{x) in (|2.8|l . Restricting (j){r,x) to the boundary will give 4'o{x). This is of 
course the standard way of doing things in AdS/CFT, and as is well-known (e.g. |21]) leads us to the 
bulk-to-boundary propagator, 



(l-/c2s2) 



• {N''^, + 2k)G'■^'■''\r,x,s,x') 



dM' 



1 



(l-fc2s2) 



K'^''-''\r,x,x') 



dM' 



(2.20) 

^"^Here we switch to the coordinates of 12. li . where the boundary is at r = 1/k, and will use the notation Xr = {r,x) and 
Xs = (s, x') for the two bulk points, where r, s are the corresponding radial coordinates, and x, x' are shorthand for the angular 
coordinates {9r,ipr,^r) and {6g,ips,tps) respectively. See section fs.ll for A'"'" in AdS and section [6. II for Af" in QTN. 

^^We found it convenient to rescale vi to half its value relative to which leads to slight changes in the formulas compared 
to that paper. 
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However we could also take a A = 2 point of view, and solve the Robin problem by directly writing a 
single-layer potential 



(^(r, x) 



dM' 



(l-fc2s2) 



•G('^''=)(r,x,s,a;') 



(j)Q{x')An{x') (2.21) 

dM' 



where here we recover (po^x) by taking {N^dr + k)(f>{y, x) and restricting to the boundary. Note that in both 
cases the same object, the bulk-to-boundary propagator (which is simply the scaled restriction of one point 
of G'(''''=)(r ) to the boundary), appears, but (j>Q{x) can be recovered from it in two different ways. Let 

us now write it down: 



s^i/fc [ (1 — fc^s^) j ^/A q^+q^—Akrcvi 

where 

q± = {l + k^r^ ± VA + 2/ifcr2v2) f 1 ± ^x/A - fi^r\2v2 - 1)) ' ^'"^^ , (2.23) 



k 
with 

A = (1 - k^rY - Ak^r^{v2 - 1) + Aii^k^r'^Y2iv2 - 1) (2.24) 

and 

c = [(1 - fi^k\^) (1 ~ mV^')] ^^'"^^ + + • (2.25) 

This Poisson-like kernel scales at the boundary as (1 — k^r^)'^S^^\x, x'y^ , which is the behaviour needed to 
obtain 00 (a;) from H2.19|l . On the other hand, {N^dr + k)K^^^■■'^'^{r,x,x') also scales as {1 — k'^r'^)'^S^^^{x,x'), 
showing that we could also recover (j>o{x) from H2.21|l ^'^. So indeed both ways of looking at this problem are 
equivalent. Of course, in both cases we obtain the same boundary behaviour, and thus the same conformal 
dimension for the dual operator. Given the scaling of 4>q{x) it should be a density on the boundary of 
conformal dimension A = 2, so (as standard in AdS/CFT) it will be a source for a dual operator O of 
conformal dimension A = 1. We will obtain the two-point function of O in the next section, and will 
eventually give an independent verification that it is conformally invariant and the dimension of O is indeed 
1 in section 121^^ 



2.3 The boundary two— point function 

We can now follow standard AdS/CFT procedures |3] to obtain the two-point function on the boundary. 
We evaluate the action (|2.4|l by first integrating by parts and then imposing the bulk equations of motion 
((V^ — l/6TZ)(j) = 0), which results in a pure boundary term: 

S=-l- lim / d''xy^{N^)<j>ir,x)(N-dr-^lc)cl,{r,x) . (2.26) 

Zr^l/kjgj^ V •J / 

Now we can substitute the solution for (j>{r, x) in terms of the bulk-to-boundary propagator 

(j){r,x)^ [ K'■^''''\r,x,x')M^')'i^ix') (2.27) 

JdM 

(where df2(a;') is the measure with respect to the finite boundary metric) which gives 

{N-dr - ^/C)0(r, x)^{l~ k^YM^) , (2.28) 

'^^By this we mean that it behaves roughly like linij^o t / {^^ + x^) —> e^S^^^ {x, x') close to the boundary. Properly normalising 
would require integrating it over the Berger sphere, which is a rather involved calculation and we have not performed it. 
^^The difference is that now one obtains t^S^''^^ (x, x') from a distribution of the type lim^ ^ + x'^)^ , rather than the 
more "exceptional" case above. 

^*Note that to solve the problem of prescribing f{x) in 12. 8i . which would lead to correlation functions of a dual A = 2 
operator, we would have to find a Green's function on QTN scaling as (1 — fc^r^)^ on the boundary (which could of course 
then be interpreted as the Dirichlet Green's function for A = 2 or the Neumann Green's function for A = 1). This could be 
constructed in a similar way to but has not yet been attempted. 
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and given the scaling of the induced measure (Vh = ^(N^) ~ (1 — fc^r^) we have 

5^/ / M^)( lim ^'^^f^^'^flf^ ^M^ym^Wx'). (2.29) 



dM JdM' \r^l/k (1 - fc2r2) 

Taking the hmit r l//c, we find that the answer is finite, which would not have been the case had we 
forgotten the extrinsic curvature term^^. Note that the way we obtained H2.29|l is exactly the opposite of 
the way it is usually done in AdS/CFT, i.e. it is the normal derivative term that gives us (pQ, while the 
restriction of (p to the boundary brings in the boundary value oi K'^^^'^^r, x, x'). As discussed, we would have 
obtained the same result working in the standard way, but at various places we will find it best to think of 
4i{x) being given by a Robin problem rather than a Dirichlet problem. 

Interpreting the boundary values 00 as sources in the dual field theory, coupled in a conformally invariant 
way to an operator O and, using the notation M(x,x') := {O{0r, 'Pr,''pr)0{9s,ips,ips)) we arrive at the 
following proposal for the two-point correlation function of the boundary operator O dual to the bulk 
conformally coupled scalar <f>: 

M{x, x') = ^ a+-a_ ^^.30) 

V(5 a+ + a- — 2coVi 

where 6 = {1 — V2)(l — fi^ /k'^\2), and 

a± = (1 ± ^v2)(1 ± 2^V5 + ^(1 - 2v2))M^-^) 

2 fc fc (2.31) 

co = (l-^)^(i--A-)(l + ^). 

This slightly complicated object is what AdS/CFT tells us should arise on the boundary as a conformal 
two-point correlation function of the scalar operator 0(x) dual to the bulk conformally coupled scalar a:). 
In particular, as discussed extensively in the preceding section, we expect 0{x) to have conformal dimension 
A = 1, and this is the statement that we would like to check in the following. Of course, in the absence so 
far of a concrete string/M theory embedding (in the sense of a decoupling limit of some brane configuration) 
of this model, it is not yet clear what the dual boundary theory actually is, and thus which fields constitute 
the gauge-invariant operator 0{x). It is expected l^j that one can think of it as some (non-supersymmetric) 
deformation of (the infrared limit of) M — % SYM on round S'^ (here we are thinking of QTNxX^ (with 
X some Einstein manifold) as a deformation of AdS4 x S''), and in particular it is expected to at least contain 
scalar operators with A = 1 (since Af = 8 SYM does). On the other hand, any CFT on S^ (regardless of 
the fundamental degrees of freedom) which contains A = 1 scalar operators could, according to AdS/CFT, 
lead to correlators with the functional form of (|2.3Q|I (in the same sense that any CFT in flat space will lead 
to powers of l/(|a; — a;'|))^^, so it is clearly important to make sure that this is what is actually required by 
conformal invariance on S^. We will start developing the necessary ideas in the next section. 



2.4 The round— sphere hmit 

Before we proceed, it is instructive to consider the limit /i = 0, which should bring us to completely familiar 
AdS/CFT territory. In this limit our correlator simplifies drastically, reducing to 

M"{x,x') = lim M{x,x') - — . (2.32) 
fi^o 1 — vi 

This is a two-point function on the round sphere, and it can be easily seen to correspond to a boundary 
operator of conformal dimension 1, by switching to the upper-half space (or Poincare) notation. The 

^^It is well known, of course, that this term is the first in a series of (covariant) regulators for the boundary theory in the case 
of scalars, derivable using holographic renormalisation I25II7I . Since we are working at a rather naive level, being interested only 
in the functional form of the two-point function rather than the precise normalisation, we will not perform a similar analysis 
here. 

^^However, unlike the case of flat space, we cannot claim that M(x,x') is the unique A = 1 correlation function. See the 
conclusions for a discussion of this point. 
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conversion can be found in appendix B of where the hyperbohc space chordal distance u is written in 
two ways:^''' 

^yy'^^^ ^ (upper-half space) 

(2.33) 

(i-fc^^:.)(i-fc^,.) (polar-type) . 

The bulk-to-bulk propagator in this case is simply which, properly rescaled and restricted to the 
boundary (which is y,y' for the upper-half space) gives the two-point correlation function for the dual 
boundary operator 

^^°(^'^')- (^Z^ (2-34) 

while in our polar case the corresponding limit (r, s — s- 1/fc) simply gives us (|2.32|l . From the upper-half-space 
picture it is obvious that the conformal dimension of the corresponding operator is A = 1 . 

As discussed in the introduction, we have another way to check this result without ever leaving the 
sphere, namely we can hit M'^{x,x') with any of the four conformal Killing vectors available. Choosing e.g. 
the conformal Killing vector (denoted ki in appendix 

+ sin(^) 
t; = 2sin - cos ^—-^ Oe + {d^ + d^) (2.35) 

2. \ Z J COS I 

we can readily check that 

i^vi'Xr) + ^v(xs)) T~ = -{<Jv{Xr) + (JviXs))—^ • (2.36) 

^ ^ ' ^ " 1 — Vi 1 — Vi 

(where av{xr) = cos(^) cos(^^^-j^) and similarly for (T„(xs)). Recalling Hl.l(l . we conclude, this time directly 
on the sphere, that M'^{x, x') corresponds to a A = 1 operator. Of course, due to the lack of conformal Killing 
vectors, this simple argument fails as soon as we leave the round sphere (by allowing ^ 7^ 0). 



3 Conformal Geometry 

In this section we give a brief overview of conformal geometry, which can roughly be thought of as an 
analogue of riemannian geometry, but where only angles are important, not lengths. Reviews of conformal 
geometry, containing references to the literature, can be found in |26l 1271 1^ . 



3.1 Preliminaries 

A conformal manifold Al is a pair {M, [g\), where M is & smooth d-dimensional manifold, and [g\ is a 
riemannian metric defined only up to scale, i.e. a conformal class of metrics. That is, rather than dealing 
with smooth, positive sections of 0T*A1, as in riemannian geometry, here we are considering positive sections 
of the ray bundle R{QT*A4) (i.e. if g is a representative metric, the fibres are t^g,t e IR+) Two metrics 
in the conformal class, say and that are related as = fi^^^t^, with fl a smooth function on the 
manifold, are equivalent in conformal geometry. 

In practice, rather than implementing this equivalence directly, it is more convenient to choose a repre- 
sentative g e [g], and examine how other quantities on the manifold, defined relative to g, transform under 
conformal (Weyl) change of the metric: 17^,^ — > fi^g^j,. To keep track of these transformations, we intro- 
duce the bundle £{11]) of conformally weighted functions (conformal densities) of weight w: Under conformal 
change of the metric, a section ip £ r{£{w)) transforms as^* 

g^,^n^g^,, =^ -> (3.1) 

^'^By "upper— half space notation" we mean the usual metric ds^ = (dy-^ -|- dx^)/(fc-^y^), where y is the radial coordinate. 

'^^Note that the conformal weight w as it is defined here corresponds to the length dimension of </>■ Often in physics it is more 
convenient to use mass dimensions, which we have been calling A, in which case A = —w. In this section we will use length 
dimensions which are more natural from a conformal geometry point of view. 
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One can also define tensors by tensoring £{w) with the appropriate bundles. For instance the conformal 
metric itself can be thought of as a section of f^i/(2), signifying that it is a two-tensor transforming with 
weight 2. We will only need to discuss scalars, so we refer the reader to the reviews cited above for more 
details. 

The prototype example of conformal geometry is the sphere S'^ — SO{d + 1)/S0((i), which is acted on 
naturally by the group of conformal diffeomorphisms SO{d + 1, 1). It plays a role similar to that of M'' 
in riemannian geometry. The high degree of symmetry leads to a beautiful interplay between differential- 
geometric and group-theoretic ideas, and helps in gaining insight which is often helpful in the curved case. 

One of the main objectives in conformal geometry is to find conformal covariants, i.e. objects that 
transform in a simple way under conformal change. These could be tensors on the manifold, like the Weyl 
tensor W^p^ which (with this choice of indices) is conformally invariant. There also exist objects which, 
unlike the Weyl tensor, are only covariant in particular dimensions, as for example the Cotton tensor^^, 

Cpi/p = ^pP^w - V^P^p , where P^^ = ^7^^^. - 2(rf^ 1) ^^") 

which is conformally covariant only in d — i, where it is an indicator of conformal flatness (i.e. it vanishes 
if the space is conformally flat). One of the most influential works in this field is of course the article by 
Fefferman and Graham in which they construct scalar conformal invariants involving (derivatives of) 
a d-dimensional riemannian metric by looking at differential invariants of an associated d + 2 dimensional 
lorentzian ambient metric, where the conformal group is linearly realised. These ideas have proved very useful 
in an AdS/CFT context, where, in the context of holographic renormalisation, these conformally invariant 
objects are used as counterterms in the d + 1 dimensional gravitational action, leading to a regulated action 
for the dual d-dimensional field theory. See |7J for a review and references. 

As in riemannian geometry, an important question that arises when given a conformally covariant object 
on the sphere (which, as indicated above, is the "fiat model" of conformal geometry) is whether it can be 
extended to a conformally covariant object in the curved case, that is on an arbitrary non-conformally fiat 
manifold equipped with a conformal metric. This will be important for us, since the Berger sphere is not 
conformally fiat (which can be easily seen from the fact that the Cotton tensor does not vanish). 



3.2 Conformally invariant differential operators 

In conformal geometry it is also of great importance to construct conformally covariant differential operators, 
i.e. differential operators that transform nicely under conformal change of the metric. More precisely (e.g. 
|28|): A differential operator 2?, defined locally in terms of g^,y, its inverse, and their partial derivatives, is 
called conformally covariant if it transforms in the following way under conformal change: 

g^.u ^ f^'.g^. =^ n-'^VQ'^ . (3.3) 

Note that V will act on everything to its right, including the factor of 17". The quantity (a, b) is often called 
the conformal biweight, or bidegree, of the operator V. To obtain a conformally invariant operator, we have 
to also specify which density bundles the conformally covariant operator acts between. If it acts between a 
section u G r(£(— a)) and a section v e r{£{—b)), i.e. we have 

Vu^v (3.4) 

this equation is invariant under conformal change of the metric. The most fundamental conformally covariant 
operator is the conformal (or Yamabe) laplacian 

^-^+41^1)^' ^'-'^ 
where 7Z is the scalar curvature. Under the change g^i, OF'g^y^ we have 

d+2 d-2 

Y ^Vl- — Yn— . (3.6) 



^^Also sometimes referred to as the Cotton-York tensor 1261 or the Weyl-Schouten tensor 1291 . 
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Thus, when acting between the density bundles £{^^) and £(— ^i^), the conformal laplacian is a conformally 
invariant operator. Note that, being a second order operator, it reduces conformal weight by two. This 
intuitive fact is a consequence of working with length dimensions rather than mass dimensions (which was 
discussed in a previous footnote). There are other, perhaps less well known but also important, conformally 
covariant differential operators, as for instance the fourth-order Paneitz operator in d = 4 (see e.g. jSIj for 
a review with references) which acts on scalars and has biweight (0,4). It is given by 

P4 = V9 (a^ + 2V^ (n'^-' - Ig^'Ti^ v.^ . (3.7) 

This operator has principal part A^, where A is the positive laplacian. On S** one can generalise this 
construction to operators whose principal parts are higher powers of the laplacian (e.g. |^ E21)) called 
"conformally invariant powers of the laplacian" but the question of whether they extend to the curved case 
does not always have a positive answer. For instance, the sixth-order operator with principal part A'^ in 
d = 4 does not have a curved space analogue In general these scalar operators Ps have biweight 

(^7^' ^2^) ^^'^ ^i^l take £{^^) £{—^^). If has weight w, then the corresponding Green's function 
M{x,x') = {(j){x)(j){x')) will transform as M{x,x') {x)^'" {x')M{x,x'), and the equation 

y/gPsM{x, x') 6'^'^'^ {x, x') (3.8) 

will be conformally invariant ii w = . Starting from the invariance of this equation under both diffeo- 
morphisms and conformal change, and assuming a nontrivial conformal Killing equation, one can now easily 
prove by adapting a discussion in (and also writing J7(a;) = e''^^\ and linearising). Note of course 
that the mass conformal dimension of (j) will be A = —w = 

The operators act on scalar densities, but there are also generalisations that naturally act on tensor 
densities with various symmetry properties, as well as on forms. See e.g. jE] |^ ES] for relatively recent 
discussions and references, as well as the reviews [271 128| . 

One important application of conformally invariant operators in physics has been the construction of 
non-local actions that generalise (to higher dimensions) Polyakov's two-dimensional action [^Hl in the sense 
that, upon variation, they generate the gravitational conformal anomaly (see e.g. |11II37| ). In fact the action 
constructed by Riegert j38| for this purpose contains a fourth-order conformally invariant operator which is 
the same as (|3.7|l . found by Paneitz at around the same time. 

In an AdS/CFT setting, conformally invariant operators of a slightly different kind played a crucial role 
in |39| . There the goal was to construct the conformally invariant equations obeyed by partially conserved 
tensors in the boundary theory, which are dual to partially massless fields in the bulk. 

So far we have been discussing differential^ i.e. local, operators. More relevant to our goal in this article 
is that one can also construct non-local conformally invariant operators |4U[ 14 1| . They will be discussed in 
due course, in section H31 after we first introduce the Dirichlet-to-Robin operator (which will turn out to 
be one of them) from a different perspective, that of inverse boundary problems. 

4 The Dirichlet— to— Robin operator 

In this section we first define the Dirichlet-to-Neumann operator, which plays a central role in dealing with 
inverse boundary value problems (i.e. the general issue of obtaining information about the interior of a 
domain through boundary measurements). We sketch a way (detailed in |321ES|) to explicitly construct 
the Dirichlet-to-Neumann operator for a given manifold with boundary. However, the object we need is 
actually the slightly more general Dirichlet-to-Robin operator, and the boundary in question actually lies 
at conformal infinity, so we will discuss some plausible guesses of how the construction might generalise to 
this case. Then, although the application we have in mind is very different from the usage of the Dirichlet- 
to-Neumann/Robin operator in the inverse boundary problems we mentioned, we give a short overview of 
the main ideas. After that we dicuss the crucial fact that the Dirichlet-to-Robin operator is conformally 
covariant^ and finally explain what all that has to do with the proposed QTN/Berger sphere correspondence 
and the question of the conformal invariance of our correlation function M{x, x'). 
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4.1 Definition 



Here, following |42[ I43| . we will give a definition of the Dirichlet-to-Neumann operator which is natural from 
the viewpoint of riemannian geometry. These references also discuss slightly different definitions which can be 
mapped to this one. So let us assume given a smooth riemannian metric g^j^ on a manifold-with-boundary 
and also some function / defined on the boundary dM. The Dirichlet problem for the laplacian A 
requires finding a function u such that: 

Au = in M {M is the interior M) , . ^ ^ 

u\dM = f ■ 

In words, we need to find a harmonic function in the bulk, that restricts to the given data on the boundary. 
Similarly, the Neumann problem asks for a function u which satisfies 

Au^O in M , . 

{N-u)\9M=h, ^^■^> 

where /i is a function on the boundary, and N is the outward unit normal. Now, if we solve the Dirichlet 
problem as posed in (|4.1|) to obtain u, it is a valid question to ask what is the Neumann data h defined by 
u through H4.2|l . So we define the Dirichlet-to-Neumann operator A as: 

Af = [N-u]\aM- (4.3) 

That is, the operator A takes Dirichlet data to Neumann data. 

One can thus think of the Dirichlet-to-Neumann operator as a composition of two operations: First one 
extends the prescribed boundary data / to the bulk in a way that satisfies (|4.1|) , and then restricts the normal 
derivative of the bulk solution to the boundary. The main point now is that A can usefully be considered 
as an operator directly on the boundary, which acts on the boundary function / and gives another boundary 
function h as the result. Note, however, that since the bulk equation admits two independent solutions, 
one of which gives the Dirichlet and the other one the Neumann data on the boundary, from a boundary 
perspective this data cannot be locally related. Thus we should expect that the Dirichlet-to-Neumann 
operator, if viewed as a boundary operator, will be non-local. Indeed there are theorems (see e.g. |44| ') 
which assure us that this is a pseudodifferential operator on the boundary, of order one, and that in fact its 
principal symbol is the square root of the boundary laplacian^*^. One (iterative) approach to obtaining its 
full symbol for a smooth riemannian manifold has been outlined by Lee and Uhlmann in |42l I43j : Start by 
writing the d + 1-dimensional metric in boundary normal coordinates {x^, ..,x'^,y}, where y > and the 
boundary is defined by y = 0. The remaining coordinates {x^, ...,x'^} are local coordinates for dA4, and 
smoothly extend to coordinates on a neighbourhood in of a point p G dA4^^. In these coordinates the 
metric is written as 

9t,i, = + {y, x)dx'dx^ (4.4) 



and the laplacian takes the form 



1 



V' dyy + -h'^dyhijdy + h- ^ O^k^ k'^ 8 j . (4.5) 

The last term clearly reduces to the boundary laplacian (which we call A^) when y = 0. Now the crucial step 
is the following j42j : Suppose we can find a pseudodifferential, order one operator A, depending smoothly 
on y, such that the Laplacian factorises in the following way: 

= {dy + E{y,x)+A{y,x,dMdv-A{y,x,d^)) (4.6) 

where E{y,x) = ^h^^dyhij. Then A, restricted to the boundary, is essentially the Dirichlet-to-Neumann 
operator^^. The proof (which we do not give here) is based on showing that {dy—A{y, x, dx))u is a smoothing 

^''See appendix |B] for some basic definitions on pseudodifferential operators. 

■^'^See also R151 for a discussion of boundary normal coordinates in a similar context. 

•^^In I42II43I . the Dirichlet-to-Neumann operator is defined as taking functions to d-forms, so the definition is slightly different. 
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operator (an operator "of order — oo") on the boundary, so, modulo smoothing operators, the relation 
dyu\dM = ^(0, a:,^^:) • f{x) holds on the boundary. See also , Ch. Ill for a related discussion. 

As mentioned, this method allows us to construct the Dirichlet-to-Neumann operator in an iterative 
manner. We can start by expanding the full symbol a of the operator A as 

oo 

(4.7) 

i=0 

where each term ai-i in this expansion is the symbol of an order-(l — i) operator. Taking the symbol on 
both sides of (|4.6|) . expanding, and equating the highest degree terms, it follows easily that ai, the principal 

symbol of A, is given by 

ai = -Vf^2(Ai3), (4.8) 

where (J2{1^b) denotes the principal symbol of the boundary laplacian, and the choice of sign corresponds 
to the outward unit normal vector giving the positive square root, 

N ■ u\qm = Kf{x) {-dy)u{y, x)\y=o = iV^ +■■■ )f{x) • (4.9) 

Here the "• • • " corresponds to the lower order terms in A. Substituting this result back into (|4.t)|) we can get 
a recursion formula for the ai_iS, and so in principle we can determine all the terms in the series (|4.7|) for 
A. There will obviously always be an ambiguity since we can arbitrarily add operators whose symbols are 
in the smoothing operators alluded to above. 

For reasons to be fully explained in section ^31 we will actually be interested in the Dirichlet-to-Robin 
operator. This can be defined by analogy with 14.3|l . as the operator taking Dirichlet data to Robin data: 

Bf^[iN + a)u]\9M (4.10) 

where a is some function defined on the boundary, and we have named the operator B to distinguish it from 
the Dirichlet-to-Neumann operator. 

Motivated by the method we just outlined for the Dirichlet-to-Neumann operator, we can try to obtain 
the Dirichlet-to-Robin operator for hyperbolic space and QTN by factorising the bulk laplacian. There are 
several extrapolations we make that seem reasonable (and certainly work) but which do not follow from the 
discussion above and perhaps lack rigorous justification. First, of course, in the asymptotically hyperbolic 
case we are dealing with, the laplacian vanishes at the boundary. Therefore we will take the (perhaps) 
pragmatic approach of first finding A at some finite radius (i.e. the "boundary" will be a hypersurface 
corresponding to a fixed value of i — the radial coordinate on the ball) and then taking it to the boundary at 
infinity via some limiting process^^. 

Second, we will treat H4.6|l as an exact factorisation. That is, we will look for an A that satisfies H4.6|l 
precisely, without the need for smoothing operators. In practice, that means that we will not proceed 
iteratively, as in 1^ . but wiU treat (HSJ as an actual differential equation (in the radial variable). (See in 
this respect the discussion in |14Q, chapter III.) 

Third, we would like to avoid boundary normal coordinates. The reason is that writing QTN in these 
coordinates will turn the metric into a series expansion in the radial coordinate (which we have been denoting 
by y here), while we find it easier to deal with the exact metric (|2.1|l and exact laplacian. Ultimately this is 
probably a matter of taste, but a more serious issue is the fact (see appendix IB)) that the complete symbol of 
a pseudodifferential operator is not coordinate-invariant. Since we would eventually like to use the operator 
we find in polar-type coordinates, we had better calculate it in polar-type coordinates. 

A final subtlety which arises is that, as we will see in section ^31 the Dirichlet-to-Robin operator will 
turn out to be a Robin derivative of a double layer potential. As is well known in potential theory, (see e.g. 
|48[ 1491 \^). although the single layer potential is continuous at the boundary, the double layer potential 
typically exhibits a discontinuity. This means that the boundary data of the double layer potential is given 
by the (appropriately normalised) difference of the limits of the solutions to the interior and the exterior 

^^Note here that the methods of |3 so far have only produced the solution of the Dirichlet problem at infinity of QTN, not 
at some finite radius. Thus, unlike the AdS case where (thanks to the ability to Fourier transform) the equivalent of M{x,x') 
is also available at some finite radius |47J , in the QTN case we know it solely at infinity. So, although presumably A. at finite 
radius contains useful information, we will need to take its limit at infinity to compare with M(x,x'). 
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problems. The Neumann derivative of a double layer has no jump, but since we are taking the Robin 
derivative, i.e. we are adding a constant term, it is not unreasonable to expect that the Dirichlet-to-Robin 
operator we are looking for will also have to exhibit this discontinuity. To cover this possibility, the final part 
of the practical recipe we give is to find A by factorising the laplacian on the inside and on the "outside" of 
QTN, i.e. the regions r < 1/k and r > 1/fc, and then taking the difference of the limits r — > 1/fc on both 
sides. 

Needless to say, the fact that we might need to examine the exterior boundary problem in order to 
properly define the Dirichlet-to-Robin operator, although very natural if we think in terms of boundary 
value problems on the ball in IR'', is rather unusual from an AdS/CFT point of view, and potentially 
raises some conceptual issues. Hyperbolic and asymptotically hyperbolic metrics on the ball blow up on the 
boundary, which we thus interpret as the "conformal infinity" of the manifold. It is thus very natural to 
discuss holography in terms of only what is happening in the interior. We will see later that this distinction 
is actually irrelevant for hyperbolic space, which possibly explains why no need for it has been seen so far 
in discussions of AdS/CFT. It could be that some new principle is at play in the asymptotically hyperbolic 
case, or (more likely perhaps) that our rough analysis above can be refined to avoid having to look at the 
exterior. For now we will just blindly accept this recipe and see that it gives the expected answer. 

Another unpleasant feature is that QTN is generically incomplete for r > 1/fc, containing a conical 
singularity. A brief analysis can be found in ySlj , where it is shown that the metric can be extended to a 
smooth one only for certain discrete values of /i. (See also [22 and [S3] for related matters.). This is probably 
an important issue, but for the time being we will simply ignore it, by assuming that since the singularity 
is at r — » 00, and we are only interested in some region close to the boundary r = 1/fc on both sides, it will 
not affect the limiting procedure. 

It should be clear that the prescription we gave above, although perhaps reasonable, is rather ad hoc. 
In this respect we should mention that there exists a formalism which has been explicitly constructed to 
deal with boundary value problems on asymptotically hyperbolic manifolds, namely the 0-pseudodifferential 
calculus of Mazzeo and Melrose [281 1541 Some recent applications (in the context of scattering theory) 

to asymptotically hyperbolic manifolds are j57ll58[IS^I6U) . Part of the idea is to treat operators that vanish 
on the boundary, e.g. the laplacian, by introducing a so-called "stretched product" which is a blow up of 
the diagonal in x (the space that Green's functions naturally inhabit) at the boundary dAi x dAi. 
We expect that working within this formalism will resolve, or at least clarify, the various problematic issues 
we cited. Regrettably, we have nothing more to say on these matters here. 

4.2 Inverse boundary value problems 

Dirichlet-to-Neumann and Dirichlet-to-Robin maps are extensively used in the discussion of inverse bound- 
ary value problems, a field that has its roots in a very influential paper of Calderon and, as mentioned, 
concerns obtaining information about the interior of some domain from suitable measurements made at its 
boundary. In particular, one could ask the question: Assuming knowledge of voltages and current fluxes at 
the boundary of an object, can we infer the (generically anisotropic) conductivity of the bulk? Although the 
main motivation for these problems (which come under the name of "Electrical impedance tomography" ) 
arose in geophysics (understanding properties of the earth's interior by measuring at the surface), and also 
in medical physics (non-invasively determining the density function of e.g. a human head by measuring 
the response to an applied current on the skin), they can be precisely mapped to questions on riemannian 
manifolds that are interesting in their own right. For more information and references, we suggest the re- 
views 2niE21 and the book • One of the main results in this field 02] is that in d > 3, and for {Ai,g) 
a compact, strongly convex, real-analytic riemannian manifold with boundary, the Dirichlet-to-Neumann 
operator uniquely determines the metric in the bulk, in the sense that if there is another metric g' such 
that Ag — Kgi , these two metrics are diffeomorphic, with the diffeomorphism equal to the identity on the 
boundary^"*. 

It goes without saying that inverse problems of this kind become especially relevant when embedded in a 
holographic context, since they are related to the question (e.g. [SSJIM]) of how information about the bulk 

^*The proof of this is based on the expansion we discussed in section l4.ll Knowledge of the coefficients ai_i is enough to 
find all terms in the Taylor expansion of the metric in boundary normal coordinates in a local neighbourhood, and (given the 
convexity requirement) this can be extended to the full metric 1421 . 
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is encoded in the boundary theory, and how to extract it ( "decode the hologram" ) from boundary data (hke 
correlation functions). The recent article j65| surveys many of the known results and to what extent they 
apply to manifolds that might be interesting from an AdS/CFT point of view. We direct the reader to that 
article for discussion and references. 

A related field that seems to be a natural setting for discussing holography is that of geometric scattering 
theory (see e.g. [3^] for a survey). In geometric scattering (and the subfield of inverse scattering) the 
main object of study is the scattering matrix 5(A). In the context of scattering on hyperbohc (e.g. {66j ) 
and asymptotically hyperbolic manifolds (see e.g. |57[ I58L 15^ l6Up . on which we recall that the solution of 
A — X{d — A) is given by 

fiy,x)y^~^ + h{y,x)y^ (4.11) 

the scattering matrix is defined as the operator taking / to h. It thus clearly generalises the Dirichlet-to- 
Robin operator, being a pseudodifferential operator of order Re(A) — d/2. Note that usually one discusses 
scattering in the "physical scattering" region, where Re(A) = d/2 but A has an imaginary part (in other 
words, in the region < — o?^/4, where the continuous spectrum of the laplacian lies), but it can be 
analytically continued to a large part of the complex A-plane, and thus describe "boundary value" type 
problems. As noted at the end of the previous section, in the context of geometric scattering there is a 
well-developed formalism for dealing with differential operators and Green's functions on asymptotically 
hyperbolic metrics. 

Returning to holography, the interpretation of the Dirichlet-to-Robin operator (and, more generally, the 
scattering matrix) in this context (but from a direct, rather than inverse, viewpoint) is easy to understand: 
In the dual gauge theory, taking Dirichlet data to Neumann/Robin data has the effect of going from sources 
to vacuum expectation values. So the action of the Dirichlet-to-Robin operator corresponds to "taking the 
path integral" in the gauge theory. Although perhaps there is no real content in the above statement, it 
is fair to say that there exist several unexplored connections between the mathematical field of scattering 
theory and the AdS/CFT correspondence (and its generalisations to more general manifolds). However in 
this article we have in mind a very different application of the Dirichlet-to-Robin operator, which we finally 
start discussing in the next subsection. 



4.3 The Dirichlet— to— Robin operator as a conformally invariant operator 

The discussion of the Dirichlet-to-Neumann operator in section 14.11 can be embedded in the framework 
of conformal geometry (which was the topic of section OJ, based on the following observation (e.g. |27j): 
Let (A4,g) be a riemannian manifold, and consider boundary value problems for the conformal laplacian 
Y (defined in (|3.5|l above). Let us take the dimension of the manifold to be n in this section. It is often 
convenient to define a boundary value problem as a pair {V, B), where V is the bulk operator under study, and 
B is the operator imposing the boundary condition. The boundary value problem is conformally covariant if 
both V and B transform covariantly under conformal change, and also if the initial weights of both operators 
match. For the conformal laplacian, there are two well-known possibilities for B that respect conformal 
covariance The Dirichlet operator 

Bd ■ u\aM (4-12) 

and the conformal Neumann, or Robin, operator^*^ 

Br: u^(^N- ■ (4-13) 

The conformal biweights of these two operators are (a, a), where a is arbitrary (restriction to the boundary 
does not affect the conformal weight of a density), and (-^^-j^, §) respectively. Since the conformal laplacian 
has biweight (^3^, ^^), both {Y,Bd) (choosing a — ^3^) and {Y,Br) are conformally covariant. 

If the Dirichlet and Robin problems for Y are conformally covariant, it follows that the Dirichlet-to- 
Robin operator that takes Dirichlet data to Robin data should also be covariant. Following Branson |223, 
we can show this explicitly by starting with Dirichlet data, given by a smooth function / in dM, such that 
its extension to the bulk Bj^^f is smooth and satisfies YBj^^f — 0. {B]j^ can be thought of as the operator 

■^^As always, we use "Robin" as short for "conformal Neumann". 



16 



taking / to its double layer potential, and is clearly unique.) Since the expression YBjj^f = is conformally 
invariant, the covariance of Y implies that B]^^f — > fl~^'^~^^^'^B]j^ f under conformal change, and thus 

B^\n-^f) = n-'^{Ba\f) (4.14) 

which means B^^ has biweight (^^^ IL^)^ which is exactly what is needed for the Dirichlet-to-Robin 
operator B — Br ■ B]j^ to be conformally invariant. It also follows that B is an operator of conformal 
biweight (^^, f ). So it acts invariantly between the bundles £dM{^ — n/2) and n/2), reducing the 

conformal weight by one, as should be the case for a first order operator. 

4.4 Implications for holography 

Let us now return to the problem of checking our proposal H2.3U|I for the boundary two-point function on 
the Berger sphere. We saw that the Dirichlet-to-Robin operator defined as above on a boundary manifold 
of dimension d — n — I acts naturally and in a conformally invariant way on boundary densities of conformal 
weight 1 — n/2. Restricting to n = 4, we see that it acts on densities of conformal weight —1 on the 
three-dimensional boundary. 

The next crucial bit of information we will need is the fact that, on the boundary, the Green's function 
of B is nothing but the (restriction to the boundary of the) Robin Green's function. We will show this by 
following, step-by-step, a very simple argument in j67| . adapting the notation in that article to suit our 
needs. Let us start by thinking of the solutions of the Dirichlet and Robin problems (for an arbitrary smooth 
manifold M with boundary dM) in terms of single and double layer potentials. In this approach we assume 
known the corresponding Green's functions GD{y,x,y' ,x') and GR{y,x,y' ,x') in the bulk^^, which behave 
as (see the discussion in section 2, and |21|) 

Goiy, X, y\ x-') ^ y^ , Gfl(y, f , y\ x') ,[N+{d- X)]GR{y, x, y' , x') ^ 0{y^) , (4.15) 

where N is the outward unit normal vector. As in section [2. 21 we write for a function u{y, x), satisfying the 
laplacian in the bulk, the expression 

uiy, x) ^ /(f)/-^ + hix;)y^ + ■■■ (4.16) 

In writing these formulas we are mainly thinking of the cases A = 1 and A = 2, where the two independent 
solutions are found at contiguous powers of y. It is convenient to introduce the boundary operators 

B^o ■■ u{y, x) ^ [y^-My. x)] and B^r : u{y, x) ^ [y~^ {N + {d - X)) u{y, x)] . (4.17) 

They are designed to either restrict to the boundary, or to restrict to the boundary after applying the Robin 
derivative, respectively. Note that we have included the correct scaling factors which will give a finite result 
acting on u as in H4.16|l . so B^ ■ u{y, x) = f{x) and likewise B'^ ■ u{y, x) ^ {d — 2X)h{x). Now let us restrict 
to d = 3, A = 2, and thus write simply Bd and Br from now on. In that case, as discussed in section [2.21 
the bulk Green's function (|2.13() corresponds to GR{y,x,y',x'). 

As usual, the Robin problem (specifying h(x) on the boundary) can be solved via a single-layer potential 

u{y,x!) = - f {BD)y'GR{y,x,y' ,x')h{x')An{x') (4.18) 

which clearly, given (|4.15(l . gives (BR)y ■ u{y,x) — ~h{x). However, we can do something different: We can 
restrict u{y,x), given by (I4.18|) . to the boundary, and thus obtain the Dirichlet data /(a?) as 

f{x) - - / {BD)y{BD)y'GR{y,x,y',x')h{x')dn{x') . (4.19) 

JdM' 

^^In this section we use coordinates (j/, x), where y is the radial variable and x = (x^, . . . , x"^) parametrise the d-dimensional 
boundary of a conformally compact manifold M. If M is AdS, the metric would be the usual "upper— half-space" metric. For 
simplicity, we set fc = 1. Note finally that the outward normal vector in these coordinates is roughly A'^ = —ydy (and exactly 
that for AdS). 
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This is a boundary operator, relating f{x) and h{x). Note that the kernel in (|4.19|) is simply the appropriately 
scaled restriction of (both points of) the Robin Green's function to the boundary, which is precisely what 
we have been calling M{x,x'). On the other hand, we can definitely find f{x) in the usual way from the 
Dirichlet function GD{y,x,y' ,x'), using a double layer potential: 

u{y,x) = - f iBR)y, ■ GD{y,x,y',x')f{x')An{x') . (4.20) 

But taking the Robin derivative of this expression for u{y, x) and restricting it to the boundary, we are 
assured (by the uniqueness of the Cauchy problem) to get back h{x): 

h{x) = -{BR)y ■ u{y,x) = f {BR)y{BR)y, ■ GD{y,x,y',x')f{x')dn{x') . (4.21) 

JdM' 

So we have found a function on the boundary giving ~h{x) from f{x). Of course, by definition, this is 
nothing but the Dirichlet-to-Robin operator B, so we can write 

{BR)y{BR)y,GD{y, X, y', x'){- ■ ■ )dn{x') = S,(- • • ) (4.22) 

dM' 

as an operator^^. But now substituting H4.19|l in H4.21|l gives 

h{x)^- f f {BR)y{BR)y,GD{y,x,y\x') ■ {BD)y'{BD)y''GR{y\x' ,y'\x'')h{x'')dn{x'^^^ 

JdM" JdM' 

[{BR)y{BR)y,GD{y,x, y',x')] M{x',x")dn{x')] h{x")dn{x") 

dM" UdM' J 

B^M{x,x")h{x")dn{x") 

dM" 

(4.23) 

in other words, 

B^M(x,x') = —6^'^\x,x') . (4.24) 

This rather general argument shows that the Dirichlet-to-Robin operator is the inverse to the restriction 
of the Robin Green's function to the boundary. This is the main result that, combined with the conformal 
covariance of the Dirichlet-to-Robin operator, will be used to check whether M(x,x') has the form that is 
required by conformal invariance on the Berger sphere. 

Let us also point out that in the hyperbolic case B can be thought of in a group-theoretical way, as an 
intertwining operator acting between representations of the conformal group S0(d+1, 1) |68[l27j . In fact this 
aspect is treated in the AdS/CFT context by Dobrev where it is shown that the operator interpolating 
between / and h (called 0o and O there) is the boundary correlator G\{x — x') ^ l/{x — x'Y^: 

h{x) - / G\{x,x')f{x')d'^x' . (4.25) 



Let us actually check this using the definition of the Dirichlet-to-Robin operator, in upper-half-space 
coordinates, where the boundary is H'': In this case, as we will show towards the end of section I^^Tl the 
Dirichlet-to-Robin operator for hyperbolic space is simply 

B = ^/A^ (4.26) 

where Aq denotes the positive boundary laplacian. We can now apply the definition of the square root of 
the laplacian, which, according to appendix IbI can be given (in IR'') through the Fourier transform: 



fix')d''x' . (4.27) 



'^See also 1451 . chapter 4, for a diflferent proof of this. 
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Now we can use the formula (see also [TU], p. 67) 



Inverting (by multiplying with e~*^^ and integrating over d'^y), and specialising to d = 3, A = 2 we find 

/i^i^^^^^""'^^^-l(djw- ^'-^'^ 

This is the same as the kernel G\ in (|4.25(l . so we have verified that the intertwining operator relating the 
two independent boundary values in the group-theoretical approach of is simply the Dirichlet-to-Robin 
operator. Note that a very similar calculation can be found in |71| . 

After these remarks, let us now specialise to QTN and conformal coupling. We now know that the 
(pseudo)differential operator, whose Green's function is our correlation function M{x,x') (|2.30|l . is the 
Dirichlet-to-Robin operator B. Given the discussion above on the conformal invariance of B, we see that if 
we can verify H4.24|l for M (x, x'), we can immediately conclude that it is a two-point function of a boundary 
operator of conformal dimension 1. In this case, this is simply a nice (and, we believe, necessary) check 
of our results, since M{x,x') should satisfy H4.24|l by construction. However perhaps cases will arise where 
one is given a boundary correlation function without knowing its origins as a restriction of a bulk Green's 
function, and then H4.24|l could provide a way of checking whether it corresponds to a A = 1 operator. 

Having discussed the background ideas we need, we will now move on to the calculations. 



5 The Round Sphere 

Before applying the methods we discussed previously to the Berger sphere, we believe it is instructive to 
do so for the simplest possible case, namely the round three-sphere, which has HH^ as infilling manifold. 
This case, which of course corresponds to the standard AdS/CFT situation on AdS4, has many simplifying 
features which allow us to demonstrate the main idea of the calculation without going into details. 



5.1 Factorisation of the hyperbolic laplacian 

So let us start by applying the ideas of section^ which boil down to trying to find a suitable factorisation 
of the hyperbolic space laplacian. We begin by recalling the polar-type metric on global IHH^: 



(l-fc2r2)2 

The conformal Laplace-Beltrami operator is 

(l-fcV)2 



2k^ 



3 + 4 

Orr H Or + ^(-Aq) 



2A:2 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



where TV'' = (1 — k'^r'^)/2 is the r-component of the contravariant unit normal vector = (iV, 0, 0, 0) to a 
surface defined by r = rg (where tq is some constant), /C is the trace of the extrinsic curvature, 



Here Aq denotes the positive laplacian on the boundary three-sphere (see also appendix 

- Aq = V|3 = dee + cot 9de + csc^ 9{d^^ - 2 cos 0d^^ + d^^) . 
The right -hand-side of (|5.2(l can be easily factorised to give 

+ 2fc2 = ( N''dr -'ilC + M ( N^'dr -lie -A 



3(l + fc2r2) 
2r 



(5.5) 
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and the pseudodifFerential first order operator A is given by 



.28 



1 — k^r 



2J2 



VAo + 1/4 . (5.6) 

r 

This can be easily checked by expanding out the right-hand side of (|5.4|l and equating with the left-hand 
side. Notice that the linear in A terms on the right-hand side 

1 . ( \- Pr^ 



-N^drA + -)CA = dr j - j v/Ao + 1/4 (5.7) 

cancel on their own. This means we could have taken A to have the opposite sign and it would still be a 
solution. So we have made a choice of sign here, which is dictated by requiring that ^ be a positive operator. 
This can be traced back to the choice of the positive square root for the outer normal vector in (|4.9I) ^^. 
Now following the arguments in section 01 we conclude that the Dirichlet-to-Robin operator at the inner 
boundary (r = Tq < 1/A;) is 



A\r=ro ■ (5.8) 



The last step is to rescale A to obtain an operator on the boundary at infinity. The most obvious way is to 
simply scale to the Dirichlet-to-Robin operator B on this confornial boundary as: 

1 



B= lim -^A = fcv/Ao + 1/4 . (5.9) 

So here we use the same defining function we used to define the boundary laplacian (see appendix 0. This 
operator is the same as the one appearing (as an example of a conformally invariant operator) in Branson's 
work |SS|0n] (where the bulk manifold of study is not hyperbolic space but IR x S'*).'^" 

Although it is not needed in the following, it is instructive to also derive the Dirichlet-to-Robin operator 
for AdS4 in the Poincare or upper-half-space parametrisation. Here the laplacian is simply (see e.g. |S]) 

Vp = y^dyy - 2ydy - j/^As (5.10) 

where here A^ — — X^i'^i simply the positive laplacian on IR'^, the boundary of AdS4 in upper-half 

space coordinates. Now the conformally coupled laplacian can be easily factorised as (we set k — 1 here) 

+ 2 = (^-ydy + 2 + y VA^) {-ydy + 1 - 2/ VA^) • (5-11) 

Note that the outward unit normal vector is —ydy, and K, — ~3. Taking the limit y ^ we find 

Bp = ^A^. (5.12) 

Note the absence of the 1/4 factor, which is a good illustration of the fact that only the principal symbol of a 
pseudodifferential operator needs to be coordinate-invariant. Lower order terms do depend on the coordinate 
system used. The form (|5.12ll of the Dirichlet-to-Robin operator was used in section [4.41 in comparing with 
the results of [5^ . 

^^As mentioned in the introduction, the factor -j could be thought of as four dimensional conformal coupling ^7^0 on the 
three-sphere, where Tlo = 3(1 - k'^r^)'^ /(2r^). 

^^The sign of A also depends on the sign in 15. 4i but that is purely conventional, and indeed here it is opposite from that in 
section f4. II Even with the opposite sign, we would still choose A as in 15.61 based on positivity. 

^"^In [68). for d = 3 the expression for B is V Ab + 1, which comes about because there the definition of the laplacian on is 
such that its eigenvalues are j{j + 2). To compare (since we are allowing half-integral I) we need to take j = 21, in which case 
(setting also k = 2 which corresponds to the unit sphere — see appendix the expressions match precisely. The reason they 
match is possibly related to that HH^ is conformal to flat space, and B is conformally invariant. 
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5.2 Eigenfunction expansion of the correlation function 

Now we need to act with B on the two-point function 

M°{x,x') ■ (5.13) 

1 — Vl 

where C is some normalisation constant. Hopefully we will obtain a delta function, demonstrating that, as 
claimed in the introduction, M'^{x,x') is the Green's function of B. The problem of course is that ;B is a 
pseudodifferential operator and it is not immediately clear (to this author!) how to act with it on a function 
on S^. We find it most convenient to expand M'^{x,x') in eigcnfunctions'^^, on which (as discussed at the 
end of appendix IB)) it is trivial to act with B since we can easily find its eigenvalues. 

So first we need to expand (|5.13|) in eigenfunctions on the round sphere. Since the round sphere is the 
symmetric space SO(4)/SO(3), there are various ways to do this. For instance we could use directly the 
results of Drummond |72j , who gives a formula for the expansion of arbitrary powers of | a; — | in terms of 
spherical harmonics'^^ . However, this expansion depends on the bi-invariance of the metric, so keeping in 
mind the equivalent calculation for the Berger sphere (where we have only left-invariance) that awaits us, 
we will choose a different route. 

To begin, we need to retrace our steps a little bit. In section we found M°(x,a;') as the limit of 
M{x, x') as fj, — 0. We could instead have taken the limits differently, that is first take /i = in K^f^-'^'^ {r, x, x'), 
and then obtain M^(x, x') as the restriction to the boundary. In that case, we write 

M°(a;,a:')= Hm C ■ ^^if ("^'^^(r, x, x') = C • lim ,„ / ^, (5.14) 

where C — 2C . Now, for reasons soon to become clear, let us decide not to take the limit r ^ 1/fc until 
later in the calculation. Then, the expansion we are looking for is: 

M\x, x') = C'-J2 a/™Yj,„(x)Y'„„(x') , (5.15) 

Imn 

where I takes all integer and half-integer values from to cxd, while m, n range from — I, —I + 1, ... I. As 
discussed in appendix^ Y,m„(x) (appropriately normahsed) equals ^/2l + le~*'"'^e~™'^P^„(cos0). To sim- 
plify matters, we can choose particular values for the coordinates of the point x', and reinstate them at 

the end. We will use 6' = ip' = i/j' = 0, which gives Y„„(0) = VWTT. Also, Vl in (|2.1t)|) simplifies to 
Vl = cos(6'/2) cos{{ip + ■0)/2). Now we are looking at the expansion 

M"(x,0) = C^^V2ma,™„Yj„„(.T) =C'-^(2Z + l)az„„e-™'^e-™'^f'L(cos0) (5.16) 

Imn Imn 

where the coefficients aimn are given by (see e.g. [75117^ ^ 

^ rS^. I' r r 1 + k^r^ - 2krl i.cos + ^„, e'-^.'-'P,;.,.(cos.)dV,d„sinW(. ^ (5^7) 

Here b ~ (— l)™^"/(167r^). First we do the tp integration, by converting to a contour integral, and end up 
with: 

2n 

a/m« =47r6 lim^y / ( ^ 'l il^ ) e*^'"-")^Pi„(cos6l)d(^sin6lde' (5.18) 




with A given by (|2.24(l with ^ — 0,6' ~ 0,tp' ~ 0, i.e. A = 1 + fc**r** — 2k^r^ coad. Performing the (trivial) ip 
integration, we get: 




2n 



aimn = b lim / —= \ — ^- — P'„„(cos 6')(5™„ sin6ld6' (5.19) 



^^This is, after all, a physics paper. 

^■^Our expression 1/(1 — vi) is really just \l\x — x'\ if we define x,x' in terms of the embedding coordinates of in as 
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This is a relatively nontrivial integral, and we can now see the value of not taking the limit r ^ 1/k 
previously. Defining h = k'^r'^ , we see that the integrand is 



1 /f + /i-x/l + /i2- 2hcos9\ 



cos ^ 



(5.20) 



2^ 



A look at, for example, [73| or [73] reveals the close similarity of (|5.2U|) to the generating function of the 
P^j„(cos0), which we reproduce here for /i < 1: 



Vl-2/icos0 + /i2 ^ Y (/-m)!(/ + m)! ""'^ ' ^ ' 

with 

_ 1 - /i cos 6* + Vl - 2/i cos 6* + /i2 
t = rr : j. ■ (5.22) 

Note that this holds for m,n (and thus I) half-integer as well as integer. Setting m = n in H5.21|l we can 
easily conclude that our integrand H5.20|l can be expanded in h as 

Substituting in the integral (j5.19|) . we have 

otunn = 8TT^b lim V S,nnh' / P;„ (cos 0)P'„„ (cos 0) . siu (5.24) 

Using the well-known orthonormality property of the P^j„(cos 6) (for any m, n, and thus also for m = n): 

f Pl,^{cose)PlJcose)sm9de^ f P:„„(z)P'_(z)dz 
Jo J -1 

(where z ~ cosO) and summing over s, we conclude 



^ <^.i (5.25) 



2/ + 1 



aimn = 167r^6 lim (5m„— — — - . (5.26) 



/i^i 2/ + 1 

We are now free to take the limit h ^ 1^^ and obtain the final formula for the ai 



167r2 2/ + 1 2/ + 1 



Qimn = 167r TTTTT = ^^n— T (5-27) 

(where we substituted 6). From this we obtain 



M°(x,0) = C ■ = C • Vj„„e~^("^+"'/')p^„(cos0) . (5.28) 

1 - vi 

Reinstating 9' ,ip' ^ijj' using (an adaptation of) the composition theorem for the P4„(cos6') (see e.g. \T,\\ . 
111.4)34 

I 

e^-^if+'^)pl^{^cos9)^ e-"'^'^-'^'h-'^^'^-'^"^Plk{cos9)Plt,{cos9') (5.29) 

k=-l 



^^Of course, one can question the validity of taking h — > 1 when I5.2ijl is really only valid for h < 1. But since orthonormality 
leaves only one term in the h-expansion for each /, namely the one proportional to h" for s = I, questions of convergence should 
not arise. Of course one can also check this by first taking /i — > 1 in I5.19i and explicitly doing the integral, say for particular 
values of l,n, using e.g. the tables for the P^„{z) given in 74 . 

^*In comparing with the formulas in [73| , note that the sign of 9' is opposite to ours, so one must use the property tl^„i~^) — 
*nm(^) given on p. 123 of that reference. 
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and reverting to the normalised eigcnfunctions Y'„„(x), we finally conclude that 

M%x, x')^C-Y, _A_yJ„„(x)y',„„(x') . (5.30) 

Imn 

It is worth stressing once more that the calculation we just went through is not the recommended way of 
expanding a function into harmonics on the round sphere! We could have arrived at this simple result with 
less effort. However the Berger sphere calculation to follow (in section will be modelled, step-by-step, 
on the previous calculation, which will thus serve as a guide. 



5.3 Checking the correlation function 

Now recall that, on the round sphere, the cigcnfimctions Yj„„(x) satisfy (see appendix IXjl 

^oYL,{x) = l{l + l)Yl^{x) . (5.31) 
We then find that the Yj„„(a;) are eigenfunctions of the operator Aq + 1/4 with eigenvalues [21 + 1)^/4 

(Ao + 1/4)y5„„(x) = B^Yl^ix) = \{2l + lfYl^{x) (5.32) 



and thus they are also eigenfunctions oi B = Aq + 1/4 with eigenvalues {21 + l)/2 



BA^{x) = ^2i±}lYl^(x) . (5.33) 
We are finally ready to act with B on A/°(a;, x') (ignoring the unknown normalisation constant): 



lf°(x, x') =B, [Y, ^YJ„„(^)yL„(x') J = E 

\/mn / Iran 



tI 



B^M\x,x') =B,\y] 7T7^Y^„(x)Y:„„(x') ) = > : ^{B,Yl^{x))Y'^^{x') 

1 ^^-^^^ 
~ ^ ^ Y„„(x)Y,„„(x ) — —-=S^' {x, X ) . 

Imn V 5 

We just rederived the well-known result that the Green's function M'^{x, x') on the boundary of hyperbolic 
space corresponds to a boundary operator O of conformal dimension Xq — I. Of course there are easier 
ways to see this (two of them have been demonstrated in section [2. 4|l but we believe it was worthwhile to go 
through all the steps before tackling the case of QTN in the next section. 



6 The Berger Sphere 

In the previous section we demonstrated that the Dirichlet-to-Robin operator is indeed the pseudodifferential 
operator whose Green's function is our two-point function M^{x,x') in the case of round S^. Now we turn 
to the Berger sphere, and try to follow the calculations in section |5l as closely as possible. Of course things 
get rather complicated, but remarkably we can obtain exact answers that generalise those for hyperbolic 
space and the round sphere. 



6.1 Factorisation of the QTN laplacian 

In this section we find the Dirichlet-to-Robin operator for the boundary of QTN. We will follow the same 
procedure as in the previous section, i.e. we will look for a factorisation of the bulk conformally coupled 
laplacian of the form: 

+ 2fc2 = (^N'^dr ~^IC + A^ (^N'^dr ~^IC-A^ . (6.1) 

^^We set k = 1 here for simplicity. 
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Here the outward unit normal vector N — N^dr and the extrinsic curvature trace /C = —VrN^ at some 
radial distance r are given by 



^.^(1^ (6.2) 



and 



As in section|0] the goal is to determine the nonlocal operator A. We proceed by expanding out the various 
terms; 

+ 2fc2 ^ {N'- fdrr + {N'-drN'- ~ /CiV'^)9^ + (2/9/C2 - l/37Va^/C) - [N-'drA - 1/3/C^) - A? . (6.4) 

Comparing with H2.3|l one sees immediately that the first two terms on the right-hand side match the 
corresponding ones in V2 + 2fc2, i.e. 



4(1 - /l2^2) 



[N^Ydrr = ^ 177^^33^ -drr , (6.5) 

and 



4r(l — /i^r^) 



The remaining terms on the left hand side of (|6.4|) are 

(l-fc2r2)2 r ^2^2(2 _ ^2^2 _ ^2^2) 



(1 — fi^r^)r^ 



^0 (l_^2fc2,4) ^1 



2k^ (6.7) 



where the positive operators Aq and Ai are defined in appendix Requiring that the left-hand side of 
(|6.4(l equals the right-hand side leads us to the following Ricatti-type equation for A: 

,2 , V, , (l-fcV)2 / ^2^2(2 _ ^2^2 _;.2^2) 



N^d.A = -A^ + -^A + ^^-^ ^Ao - J + X(.) (6.8) 

where the zeroth-order part X(r) is given by 

Xir) . 2/9/C^ - l/3iV^a./C - 2k^ ^ -^^^^^f^^^^^ ~ , (6.9) 
TZr being the scalar curvature of the induced metric at the (inner) boundary, which is (see appendix^) 

T^r = 77%^^^^?^ (3 - 8a^'^' + M'fc'r^ + Vr^) . (6.10) 



The nonlinear, first-order equation H6.8|l can be solved for A in the usual manner, by transforming it into 
a linear, second-order equation. We will outline the main steps: First, we bring it to the standard Ricatti 
form: 

drA^ aA^ + bA + c (6.11) 

with 



(6.12) 



,4 



(1 - fc2r2) y 1 - ^Pk'^r 

3 + 3fc2r2 - 4^2^2 _ 7^2;,2^4 _ ^2^,4^6 6^4;.2^6 
" 3r(l - /c2r2)(l _ ^2^2)(1 _ ^2p^4) ( ■ ) 
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and 



4An 



9 - 24/i2r2 - G^^k^r^ + 8^*r* + 32n'^k\^ - T/i^feV^ - U^i^k^^ 

+ 9^2(1 _ ^2^2)2(1 _ ^2^2^4)2 

The next step is to choose A — —-^ ^"^(r? ' '^hich gives us the equation 



(6.14) 



dr 



drrw{r) — h b I drw{r) + acw{r) = . 



(6.15) 



In order to solve this, we must make a convenient rescahng of 'w{r), to reduce (|6.15|l to a known second-order 
equation. It turns out that a particularily nice equation results if we rescale w{r) = f{r)v{r), where^^ 



/(r) = (l-^2^2)i(^_^2^2^4)t-# 



1 + /i/cr^ 



2\ |\/Ar 



(6.16) 



The resulting equation for v{r) is much simpler: 
'k + 3//2fc3r4 - 4^(1 - fc2r2)VA7 



drrV{r) — 



kr{l - ^i^k^r^) 



T (\ 4(PAo + ^fcVA7-^2Ai) ^ ^ 

9Mr) vir) = . (6.17) 



The final step is to perform a change of variables: r ^ u ~ ^t^ikr'^ (where we temporarily assume ^ is 
nonzero), after which we have (note that u{l — u) — —(1 — fi^k^r*)/ {Ajs^k'^r'^)) 



u{l-u)duuviu)+ 1 - (1 + 7^)v/Ai - 2 (1 - f v^Ai) w duv{u) + 



k 



fc2 



v{u) = 0. (6.18) 



Here we have omitted an overall factor of 4/[(l — ii^k'^r'*')r'^]. We can recognise (I6.18|) as the hypergeometric 
equation 

(u(l - u)duu + [7 - (1 + a + P)u]du - a/3) v{u) = 0, (6.19) 

with 



1 



a=-- v/Ao + 1/4- |\/a7 



/3-- + v/Ao + l/4-|yA7 



(6.20) 



7=1-1 



whose solution is generically given by (e.g. \7b\ ) 

v{u) = ciF{a, (3; 7; u) + C2U^~'^F{1 — 7 + a, 1 — 7 + /?; 2 — 7; u) 



(6.21) 



However, before we proceed to solve (|6.18|) . we have to backtrack a little and consider the range of validity 
of H6.21|l : It is really valid when \u\ < 1, and that is clearly not always the case, as can be seen from the 
definition of u above. So we will need to perform a fractional linear transformation to bring it into the 
required domain. We will distinguish two cases, the interior (r < 1/fc) and the exterior (r > 1/fc) of QTN, 
which will be dealt with separately. 

^^The MAPLE algebraic manipulation program was instrumental in finding this expression. 
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The interior of QTN Here r < 1/fc, and since also /i < fc it is clear that u = (1 + /ifcr^)/(2/iA;r^) > 1, 
with equality for r — 1/k and only if /i = fc also. So the relevant coordinate transformation (7b\ is u = 1/x. 
In terms of this x — 2/ifcr-^/(l + ^kr^), and after rescaling v{x) = x°'v{x), equation (|6.18|l becomes 



x{l — x)dxxv{x) + (1 — /3 + Q!+(7 — 2 — 2a)x) dxv{x) — a{l + a — j)v{x) . 



(6.22) 



This new equation has a' — a, (3' = 1 + a — j and 7' = 1 — p + a. There is a slight subtlety here: Recall 
that we should really be thinking of our operator as acting on eigenf unctions, in which case Aq is replaced 
by l{l + 1), and Ai by n^. But then 7' = 1 — (2Z + 1) = —21, i.e. it is a (negative) integer and the usual 
solution H6.21|l does not hold (i.e. does not lead to two independent solutions). Looking at, for instance, |7S) 
we can find the correct solution for this degenerate case: 



vix) = cix^'"' F(a' - 7' + 1, /3' - 7' + 1; 2 - 7'; x) + C2F{a' , /?'; 1 + a' + - 7'; 1 - a;) 

Replacing x with the original variable r, we thus have 

J2 \ 0-a 



v{x) — Ci 



2/xfcH 
1 + /ifcf 



(notice that here we use the original a and /3 from H6.20|l ). where 



Fi(r) - ^ Q + VAo + 1/4 + I + VAo + 1/4 ~ ^v^; 1 + 2^Ao + 1/4; 



2fikr 



1 + jikr^ 



and 



(r) = F f i - v/Ao + 1/4 + VAI, i - ^Ao + 1/4 - ^ v^; 1 + (1 - ^) TaT; 



1 + /ifcf 



(6.23) 



(6.24) 



(6.25) 



(6.26) 



Putting everything together (not forgetting the extra factor of from v{x) — x"v{x)), the solution for 
is 



A' 



1 — ii'^k'^r'^ ^ 


fir) 




) Fi(r)+c- 


f 2fj.kr^ \ 


"^2(r) 


)] 


f{r){ 


f 2fj.kr^ \ 
\l+tj.kr^ ) 


Fi(r) + c- ( 


2fikr^ \ 
l+fikr^ ) 


"P2{r)] 





(6.27) 



where the constant c = cxjci. To fix this constant, we can check the behaviour for r <C 1/fc, l//i, where the 
QTN metric resembles flat space: 



As' 



(6.28) 



By factorising the corresponding laplacian we conclude that A l/ri/Ao+T/l, for small r, where the 
important piece of information is the positive sign of A. Looking at the two extreme values of c in l|6.27|l . 
we can check that 

VAo + 1/4 



A{c^O) 



VAq + 1/4 



A{c = 00) ^ -- 



(6.29) 



r r 
Thus it is clear that taking c — gives the full answer. So we will take this as our solution for the interior 
of QTN. Exphcitly, 



A 



(1-fcV^) / l~/i2fc2r 

2 V 1-/1^2 



■dr In 



fir) 



2fikr 
1 + /ifcr^ 



2 X i + \/^o + l/4-f x/zir 



Fi(r) 



(6.30) 
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The exterior of QTN We turn to the ease r > 1/fc. First note that in this case the factorisation equation 
(|6.1I) should be modified slightly, since by convention we are defining A with respect to the outward pointing 
normal vector. If we are solving the exterior problem, the sign of the normal vector will change, and so 
will the sign of the extrinsic curvature of the boundary (seen as the sphere embedded in the 4-dimensional 
manifold that is the "whole" of QTN). So essentially we need to solve (|6.1ll with the opposite sign in front 
of A. This clearly changes nothing in the calculations, we can effectively simply flip the sign of A"^' relative 

toA\ 

On the other hand, we should check whether the solution of (|6.18|) changes if r > 1/k. This condition 
translates to u > 1/2 and we should transform as u = 1/(1 — y), so that y — {u—l)/u < 1. The limiting case 
y — I occurs when ii — (opposite to the interior problem). Substituting in (|6.18|) . and this time rescaling 
v{y) = (1 — y)°'v{y), we end up with 



2/(1 - y)dyyv{y) + (l + /? + a- 7-(2-7 + 2a)y) dyi{y) - a(l + a - 7)5(y) = 



This equation is also degenerate, in a similar way to 



In this case a" = a, f3" = 1 



(6.31) 

7, and 7" = 



1+/3+Q;— 7, so 7" is not an integer, but /?" is a negative integer, being equal to ^ — y'Ao + l/4+-\/Ai' = —l- 



acting on eigenfunctions. (We take n always less than I, see appendix IXI) . Noticing that also 7" 
is integer, we can find the two independent solutions by looking at 7(3 . So we solve (|6.31|l as 



v{y) = Fair) + 



1 + /ifcr 
1 — iikr 



2x (i-f)V3T 



1 + i^ikr^ 



2 X 2^Ao + l/4 



F,{r) 



- a" - f3" 



(6.32) 



Here a is as in 



and the two hypergeometric functions are given by 



1 



Fsir) =Fi-- VA0 + I/4+ v^, - - VAo + l/4-^v^; 1 + (1 - ^)v^; - 



and 



1 



1 



^^4(0 = F ( ^ + v/A^TI7i- TaT, ^ + v/A^TT7i+ ^ VaT; 1 - (1 - ^)^/a7; ^ 




(6.33) 



(6.34) 



The important thing to note is that ^3(7') is actually the same as the second solution F2{r) of the interior 
problem, which was discarded there. We can identify i^3(r) as the correct one to take for the exterior problem 
by the following reasoning: The transition between u > 1 and u < 1 is at the line r = l/(fik) which always 
lies outside the boundary (apart from the point fi = k). So for 1/A: < ?■ < l/{fik) we could still use one of the 
solutions to the interior problem. Assuming that the only discontinuity where the solution could change is at 
the boundary r = 1/k, clearly the solution for 1/k < r < l/(/ifc) should be the same as that for r > l/(/ifc), 
and (since Fi{r) presumably should not be used for r > l/(/ifc)) that leaves us with the solution which is 
common in both regions, that is F^(r) = F2(r). Notice that this is actually a polynomial, which is in fact 
the reason it shows up in both cases. So we write for the solution in the exterior: 




dr In 



fir) 



2fikr^ 
1 + ^kr"^ 



VAo + l/4--'i\/ST 



F:^{r) 



(6.35) 



Note the extra minus sign, which comes from the definition of the exterior problem. 

As discussed in section 14.11 we can now find the Dirichlet-to-Robin operator on the boundary of QTN 
by taking the sum of the boundary limits of A^ and A'^ . 



B 



lim — 



1 



1/k 2 (1 - fc2r2) 



[A' + A"] 



(6.36) 



where the limit is understood as coming from lower values for A^ and higher values for A"^ . This expression 
corresponds to the difference of the inner and outer normal derivatives at the boundary, and the reason we 
get a plus sign is of course that we decided to flip the sign of , so that it matches the outward normal 
vector of the exterior. We will give the final expression for B in section 16.31 after we have discussed the 
eigenfunction expansion of M{x,x'). Now we turn to some interesting special cases. 
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Special Cases One special case to consider is the situation where the operators and act on functions 
that have no ip dependence. Then we can effectively set \/ Ai ~ \f~^d^, = everywhere, and we obtain'^^ 



An 



•^0 



1 — /i^r^ 



-dr In 



/oWQ 



1 — /i^r^ 



dr In 



fo{r)P^ 



1 



(6.37) 



(6.38) 



where /o = (1 — /i^r^)s (1 — fi^k'^r'^)^ jr. Of course we could also retrieve this result by working directly with 
the hypergeometric equation H6.18|l and demonstrating that it reduces to the Legendre equation when Ai is 
ignored. 

Another special case arises when we consider the /i = limit, which should make contact with the results 
in section Ism Taking the limit /i ^ of H6.30|l and (|6.35() we find 



(1 - fcV^ 



2r 



(6.39) 



Thus we have a formula for A that holds everywhere, and precisely matches what we found for hyperbolic 
space (equation l|5.6|) ). This shows why in that case we could get away with being naive and not paying 
attention to the difference between the interior and exterior problem. Note that the opposite sign of A^ was 
crucial in obtaining this result. 



6.2 Eigenfunction expansion of the correlation function 

In the previous section we obtained an expression for the Dirichlet-to-Robin operator, with which we have 
to act on our Green's function. Now we face the task of expanding the two-point function M(x, x') in (|2.30(l 
into eigenfunctions on the squashed sphere. As discussed in appendix ^ separation of variables works in 
exactly the same way as on the round sphere, which was the topic of section l5.2l Let us then write formally: 

M{x, x') = ^l-^Ani^)Y\,^{x') (6.40) 

Imn 

where Y;„„ (x) = V2r+Te-'"'^e-'"'^P^„(cos6>) for x = {9, ip, V'}, and similarly for Y^„(x'). Recall that the 
factor of \/2l + 1 is needed for the Yj„„(a;) to be orthonormal, and also that 1 = 0,^,1,..., while m and n 
range from + We can now proceed in exactly the same way as in section by picking special 

values for the coordinates of the point x', 0' = 0, = Q^ip' = 0, which can be restored at the end. Then 
Y^„(0) = ^/2l + 1. We can now expand M{9, ip, ip) in the Yj,j„(6', (p, ip) as follows: 

M(0,^,V') = ^V2rrTa,™„Y^„J0,(p,V') =^(2; + l)ai™„e-*™'^e-"'^Pi„(cosf?) (6.41) 



Ivan 



where the coefficients aimn are given by^^ 
air. 



A//(6',(p, V)e^'"'^e"'''Pi„(cos6l)dV'd^sine'd6l . (6.42) 



We thus need to calculate a triple integral to find the aimn- As in section [5.21 it will prove convenient not 
to take the limit r ^ 1/k until the very end of the calculation. That means that we will actually replace 

^^Here we use the following formulas from )75l . p. 200, relating the hypergeometric function to the Legendre functions of the 
first and second kind : 

P^(,)=(l±iyF(-l,-l;l;^], and Q;(^) = -^^I|±iL(^ + f 1 + Z, 1 + Z; 2 + 2Z; ■ ^ 



z + lj' 2'+ir(Z + 3/2)' ' V ' ' 'l + z 



where in both cases z 

3*The factor b is (-l)'"-"/(167r2), as in section ICT 
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M{x,x') in (|6.4U|) by K{r,x,x') = K^^^'^\r,x,x') from I2.22|l and then evaluate it on the boundary at the 
end of the calculation. Let us first do the integration over -0. Converting to a contour integral by defining 
w = e*'^/^, we find 

lim / K{r,e,Lp,^)e''^'t'dil^= lim { ^~ . (6.43) 

r^i/kJo r^i/fe VA V2fcr-ccos(6'/2)eW2y ^ ^ 

where we recall from section that A = (1 — k'^r^Y ~ 4fc^r^(v2 — 1) + 4/i^/c^r'*V2(v2 — 1), (?_ is given by 
g- = (l + fcV-VA + 2/iA:r2v2)fl-^%/A-/iV(2v2-l)V^ , (6.44) 



and 

c = [(1 - Ai'fcV^) (1 - mV^')] ^^'"^^ V(l + Ai/fc)(l + M^^') • (6.45) 
Notice that with our choice of 9',(p' ,tp' = 0, we have V2 = ^(1 + cos 9), so (defining z — cos 9) we get 

A = 1 + - 2{k\^)z ~ ^(fcV)(l - z2) (6.46) 

^^^'^ 1 M 

= (l + fcV-x/A + ^fcV(l-z)) (^1-^A/A-^fcVzy ' . (6.47) 

Since /i, fc, and r (which is now viewed just as a regularising parameter) appear only in the combinations 
IJ,/k and fc^r^, it is convenient to rename them as fj,/k n, and k^r^ h (as in l5.2|) . After these cosmetic 
changes, we can write 



ai 

h^l 



T / \ 2n 

7r /'Ztt 



Jo 



where now 

A = l + h^ -2hz- fi'^h^il- z^) , (6.49) 
g_ = f 1 + - \/A + /i/i(l - z)) f 1 - ^\/A - /i^/iz) ' (6.50) 



and 

C= ((l-Ai2/i2)(i_^2))i(l-M)y(YT7I)(TT7^- (6.51) 

Note that the range of /i now is from to 1. Noticing that ip appears only in the exponential, it is trivial to 
do the integration over ip: 



We are now left with the nontrivial integral over z. Although there may be much more elegant ways of 
calculating it, motivated by the resemblance of the integrand to the generating function for the P,j„(2) 
(j5.2,S|l we chose to expand in the parameter h and write each factor in the expansion in terms of the P^„(z), 
so that we can use their orthonormality properties to perform the integration. This is the tactic we employed 
in section 15.21 but recall that there the integrand was exactly the same as the generating function. Here 
things will turn out to be somewhat more complicated. 
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A simple case: n=0 Let us first look at the case where rt = 0, in other words we want to calculate the 
coefficient ajoo for any I, where of course now / has to be integer (since n is). If we set n — 0, the integrand 
(apart from the P^ni^) factor and the Sn'^b) is just (from H6.49|l ') 

/ := ^ = (6.53) 

^l + h^-2hz-fi^h^{l-z^) 

and we expand in h as follows: 

oo 

I = ^h'Gs{z). (6.54) 

s=0 

To be concrete, let us write down the first few Gs(z)s: 

Go = 1 ; Gi = z ; G2 = -^(1 - + A^'^' - Sz^) ; G3 = -|(3 - 5z^ - 3^' + 3^l^z^) . (6.55) 

We would like to express the Gs{z)s in terms of eigenfunctions of the boundary laplacian, i.e. functions 
Hs{z) that satisfy (see appendix 

/\bHs{z)^ s{s + l)H,{z) (6.56) 

(since n = for the moment, the term proportional to 11? coming from the extra part in As will not 
play a role. That is, for now Ab=Ao.) Clearly the Hs{z) will have to be proportional to the Legendre 
functions Ps(z), and we write Hs{z) = NsPs{z), with Ng a normalisation factor to be determined. First we 
find linear combinations of the Gs that are eigenfunctions, and call these combinations Hs{z): 

Ho{z) = Go{z) ; Hi{z) = Gi(z) ; H2{z) = G^iz) - ^fi^Go{z) ; Hsiz) = G^iz) - ^/i'Gi(z) 

6 3 in 5 (6.57) 

H^iz) = Giiz) - -fi'G2{z) + ^Go(z) ; H^iz) = G^iz) - -li^G^iz) + ^fi'Gi{z) ; . . . 

and so on. It is now easy to find expressions for the Gs in terms of the Hs at each (fixed, low) order in h by 
inverting these relations. Here are the first few: 

1 3 

Go ~ Hf) ; Gi — Hi ; G2 = H2 + ; G3 = H^ + Hi 

6 1 10 3 (^-^^^ 
Gi^Hi + -fi^H2 + -Ho ;G5^H5 + + ^m'^^i ; • ■ ■ 

7 5 9 7 

As discussed, the Hs{z) are completely specified apart from their normalisations Ng. Let us write the first 
few of those: 



(6.59) 



We now have to start summing up these series to end up with a closed form. By examining the terms that 
occur for higher s, we eventually find that the series for Ng has a simple generating function: 

As for the Gs{z), it turns out that it is most convenient to treat the s even and odd cases slightly differently in 
the beginning. We eventually find (this is probably the most nontrivial part of the paper from a calculational 
point of view, but as there is little insight to be gained from the intermediate steps, we spare the reader the 
details.) 



G2r{z) = Tlf^l^'H2r^2^{z) (6.61) 

i=0 
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and 

r-l 

G2r-l{z) = ^r,(°Vj^2r-2.-l(^) (6.62) 

where the coefficients and T° of the even and odd terms respectively are given by 

rl 



(6.63) 

(6.64) 



^ l (4r + l-4i)!! (2r-l)!! r\ 
~ i! (4r + 1 - 2\)\\ (2r - 1 - 2i)!! (r - i)\ 

and 

(„) ^ 1 (4r - 1 - 4z)!! (2r - 1)!! r! 

i! (4r - 1 - 2i)!! (2r - 1 - 2i)\\ {r - i - 1)! 
So we have succeded in writing the integrand as an expansion in the Legendre functions: 

oo / ^ r — 1 \ 

^^Y.^' [Y. T!:f^i^''N2r-2^P2r-2^{^) + ^ T^^- V^2r-2»-l/'2r-2»-l (2) . (6.65) 

s=o \j:=o 4=0 / 

Now we can proceed by first doing the integral over z in H6.52|l . We get 

f I ■ Pl{z)dz ^Y.^^^'il'^rf ^^^"N2r-2^ ^ P2r-2^{z)Pl{z)dz + 
J-1 ^ -^-1 

+ E E T;^V''A^2.-2.-1 P2.-2.-l(^)PK^)dz = 



(2r-l)!!r! \ / 2-nP,(l/^) 



{2r + l + l)!!(r - 1/2)1 J \ {1/2)1(1 - 1)!! 



(6.66) 



(2r-l)!!(r-l)! \ / 2 • /!P,(1/m) 



i2r + l)\\{r- {l + l)/2y.J V ((^ - l)/2)!n! 
Here we first used the orthonormality property of the Legendre functions, 

'^P.(z)P,(z)dz = ^ (6.67) 

which gave 2r ~ 2i — I and 2r — 2i — 1 = / in the first and second sums, respectively. Summing over the 
delta functions then gives the final result. Clearly the first sum in (|6.66(l will be there only for I even, and 
the second one only for I odd, otherwise they should be thought of as being zero. Now we have to put these 
two sums H6.66II together again. 

Manipulating the series expansion, using relations such as (2rt)!! = 2"ri! (for which we first need to 
convert the double factorials of odd numbers to ones of even numbers, through (2n + 1)!! — {2n + l)!/(2ri)!!) 
and T{2n) = 22"-70Fr(7i)r(n + 1/2) (e.g.^^), and substituting r = k + l/2 for I even, we find for the even 
I series in H6.f)6|l the result 

/ 1 \ 2 / 0Fr(; + i)M'+^ \ ^ r(fc + 1/2 + i)r(fc + 1/2 + i/2)r{i + 3/2) 
' M V r(; + 3/2)2'+i y ^ r(//2 + i)r(//2 + i/2)r(fc + / + 3/2) ^ • ^ ' 

Notice that now the range of the summation variable k is from zero to infinity. Similar manipulations of the 
odd / part of (|6.t)6|l . substituting this time r — k + {I + l)/2, give the same result. This means that (|6.t)8|) 
is actually the full solution to our integral. What's more, the series over k turns out to actually be a rather 
simple function: Noticing the similarity to the hypergeometric generating function (which we have tried to 
exhibit by writing the series in this way) and recalling that the Legendre function of the second kind can be 
expressed as (e.g ^^) 

^'(^) = faTwW^'^(2+^'2 + 2^^+2=^j ^'-''^ 
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for |2:| > 1 (i.e. fi < 1 in which certainly is the case here), we have shown that 



Km 



h^lj^l ^1 + - 2hz - fi^h'^il - z"^) 



P;(^)d^ = E^^'Q)Q'Q)- (6-70) 



This final formula is refreshingly simple, so much so in fact that it cries out for a smarter (and less onerous) 
way to obtain it'^^. However recall that the purpose of doing things the way we did was to prepare us for 
the more demanding calculation for n =/= 0. As we will see, much of the work for that case has already been 
done. 

To summarise, so far we have found one coefficient in the expansion (|6.4U|) of AI{x,x') in spherical 
harmonics H6.41|l for each (integer) I: 



aioo 



(6.71) 



The general case Now let n be arbitrary. Let us again expand the integrand of H6.52() as a series in h: 

/ \ 2n ^ 



1 



q- 



A \2c^/h^{l + z)/2 



Y.h^G^^X-{z) 



(6.72) 



s=0 



The reason for the notation G*+"(2;) is that, as we will soon see, these functions are related in a simple way 
to the Wigner functions P^^"(2;). The first few terms in the series are''^" 



0: GlJz) = \{l + ^,Y+^{l-^,)^'^^]^ ■{^-±^^ 



s = 1 : Gl+'{z) = -{n^i'z - nz - z + n) ■ GH^iz) 



(6.73) 



s = 2 : G';^+^{z) ^ -(2/^2^2 - 2(1^ z^ - Infi^z'^ - An^fi^z^ + Qz^ + Inz^ + 2r?z^+ 

,2,,2„ I o. — 2„ ci — /i~„2 I o, — 2 I o,,2 o i o„2 



+ in^fi^z + 2n^^z - 6zn - Azn^ + 3n^^ + 2fi^ -2 + 2n^ - n) ■ G"^{z) 



The expressions rapidly become awkward, so we have only written out the first three terms. We can easily 
verify that G'^„{z) satisfies the SU(2) differential equation, with m — n, and eigenvalue —n{n + 1)"*^: 



(1 - z^)d. 



2zd, 



2n^ 

l + z 



Gnni^) — ~- ^^nn(^) 



(6.74) 



So G^„{z) has to be proportional to P,"„(z). On the other hand, we find that G"+^(z) is not an eigenfunction 
of IIA.28P with m = n. As in the previous calculation for n = 0, we should consider linear combinations of 
the G"s that are eigenf unctions. In this case we define 



h: 



n+l 



(z) = G\\'^ {z) H —Gnn{z) 
71+1 



which turns out to be an eigenfunction, with eigenvalue —(71 + l){n + 2): 



(1 - z^)d. 



2n' 

l + z 



1 



(6.75) 



(6.76) 



Proceeding in this way, we can find linear combinations of the G'5^^'(z) that are (unnormalised) eigenfunctions 
of the SU(2) differential equation, with eigenvalues [n + l){n + I + 1) , and thus have to be proportional to 

is one, of course: We can first take the limit h — > 1 in / and then integrate it with a Legendre function, which can 
be done for particular values of I since the Pi{z)s are polynomials. Then we can try to find a general formula for the result, 
which again leads us to l6.7Ui . 

^'^Note that this expansion holds for both integer and half-integer values of n. 

^'^In 16.741 we have added the extra factor of 1/4 which makes the eigenvalues complete squares. 
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the P^^^{z). Let us write down the first few: 

^nn (^) = ^nn (^) + ^ _|_ -j^ ^rm(^) 

^nn l^J-Wm + „ + 2 ^ ' 2(n+2)(n + 3) ^nnl^J • 

Tfie relations rapidly become unwieldy as I increases, so we will not write any more of them out explicitly. 
In any case, what we are more interested in is the inverse relation: How to write the G"s in terms of the 
i?'s. Let us again write down only a few terms: 

rn+2,._ Trn+2,. x , 2 2^ V + Zti^ + n V + n V + + 5n + 2 

^ ^~ ^ ' n + 2 ^ '^^ 2(2n + 3)(n+l) ""^ ^ (6.78) 



n + 3 ^ ' ' ^ 2(n + 2)(2n + 5) 

fi^n'^ {26 + 42n + 25^^ + 4/^2„2 g^3 _^ 3^2^3 ^ 2^2n'*) 



6(n+l)(n + 2)(2n + 3) 



We have thus managed to write the first few terms in the expansion in terms of the H^^^{z)^ which are 
just unnormalised versions of the P^^'^{z). The normalisation factors can be obtained by dividing them, 
^ snn=H'!^^^ {z) / P^n'^ [z) (whcrc wc can find the P^^^{z) either by looking at tables, or by expanding the 
generating function H5.21|l for m ~ n). We have 



_l + n — jjL n 

J^lnn — — ; • jVonn 

n + 1 

_ 6 + 7n + 2n2 - 2^"^ - l^i^n - A^i^n^ + 2^'^n'^ 
[n + 2)(2n + 3) 

_ 30 + 37n + 15^2 + 2n^ - ^^(18 + 45n + SOn^ + 6^3) + ^'^{8n + ISn^ + 6n^) - 2iJ^n^ 

(n + 2)(n + 6)[2n + 5) 

(6.79) 

Let us now write the integral 1)6.52(1 as 



f I-PL,{z)Az^Y.^^ f Glt^{z)Pl^^{z)dz^Y.^'^Nl^^Y.a^\^i) f (z)Pi„(z)d 



(6.80) 



where a*' (/it) symbolises the various coefficients in ((6.78|l . This integral receives contributions from all 
s > l — n, since then for i = n + s — lwe have a nonzero contribution. Now notice in H6.78|l that the terms to 
the right (higher i in H6.80|l 'l have higher powers of /i^, and that the power increases by fi^ every two terms 
or so (this is really only evident in G^^^{z), but can be checked to hold for the higher terms). So, truncating 
the sum in s should give us the answer to some order in /i. For example, we can check the case / = n + 1, 
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to order /i'*:'*^ 



1 



= A'i„„— , ^ 1 — - + ;T7 T^^TiT. ^ (D.Sij 



'2(n+l) + lV " + 2 2(n + 2)(2n + 5) 

2 2(2n-'' + 8n2 + 9n + 6) 



2n + 3 (2n + 5)(2n + 3)(n + 2)(n+ 1) 



So this is a contribution to when / = n + 1 . It can be easily checked that it reduces to the aioo above 
when n = 0, as is the case for ah the previous formulas. Now the question is whether we can extend these 
relations to arbitrary /, i.e. write the in closed form. This is where the special case we considered 
earlier comes in especially handy. Recall that we found that a;oo — I I^Pii^ I lAQii^ / ■ Clearly the formula 
for n 7^ must be a generalisation of this. We can cheat a bit by looking at the expressions in section IHTI 
where, as we saw, there appeared some hypergeometric functions that reduced to Legendre functions in a 
special case. Motivated by this, we can make a guess as to how to generalise Qi{\/ii). The guess 

is (see a previous footnote): 

Pi i-^ ^ f^TT^l ^ f(-1 + n, -l-fin;l + {l-fi);^'^ 



, u I \ 2ii I \ 1 + n , , ^ 

Qi (-\ ^ f f(i + 1 + 71,1 + 1~ /in; 21 + 2; 

Trying this out, and fixing coefhcients to match with the lower order terms for the aimn (that we can 
calculate as in H6.81|l l we end up with: 



, • 1 - V^'^^" r(l + ^ + n)r(l + I - fin) 



2 n(l - /x2) V 1 + r(2/ + 2)r((l - fi)n) 

X F [-1 + n,-l- fin;l + n- fin; ^ -] x F ( 1 + I + n,l + I ~ fin;2l + 2; '^^ 



(6.83) 



l + fij V r , . ^^^^ 

where we have also substituted the factor b from (|6.52() . This proposal can be checked to give the correct 
answer for the coefficients aimn that we can calculate as in (j6.81(l by expanding in fi. Note that despite 
the explicit appearance of n to the first power, its expansion in fi contains only powers of n^, as it should 
by looking at e.g. H6.78|l . Although one could probably perform further checks of this result, we take the 
matching with the low-order terms as sufficient evidence to take it as our final answer 

We can now proceed (as in 15.211 ') to restore the angles 9', ip',tp' that we had set to zero, and conclude 
that the expansion of M(a;, a;') is given by (|6.40|l . with the aimn given by H6.83|l without the 6mn- 

6.3 Checking the correlation function 

In section 16.11 wc found two different solutions for the Dirichlet-to-Robin operator, one that should hold 
"inside" the manifold (i.e. for r < 1/k) and another that should be used "outside". We also argued that the 
correct operator to consider on the boundary r = is the sum of these two: 

A^^{A'+A^). (6.84) 

Here we finally check the results of 16.11 against the eigenfunction expansion for M{x,x') that we obtained 
in the previous section. First we focus on the special case where A is acting on a function that does not 
depend on the angle tp. That means we can effectively set V^i = everywhere, and we obtain (see the end 
of section IHj]! 




/o(r)Q 



(6.85) 



Taking also /i — > 1, as in all expressions that follow. 
*^To check that we recover 16.711 for n = 0, note that nr((l — lJ,)'n) 1/(1 — /i) for n - 
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and 



1 — /i"^r^ 



dr In 



(6.86) 



where /o = (1 — yrr ) & (1 — fj. k r jr. Notice that although these are supposed to be operators on the 



boundary, where r = 1/k, we have kept r free since it will prove useful in the calculations below' 
(j6.84|l we obtain 



,44 



(l-fc^r^) /l-/i2fc2j.4 
4 V l-/iV2 



9r In 



p 



From this expression, using the property dz h\{x/y) = {ydzX — xdzu)/ (xy), and the identity^^ 

1 



Pi (z) dzQi [z) - Qi (z) dzPi {z) 



1-Z2 



we finally obtain 



A = 



(1 — k^r^) ^kr 



(1-^2^2) 2 



From 



(6.87) 



(6.88) 



(6.89) 



Now we can take r — > \/k, with the usual defining function 

1 . 



lim 



r^l/fe (1 - fc2r2) 2fc(l-^2/fc2)^ 



^A„ + l/4-i *5v^o + l/4-i 



(6.90) 



Recall that ^("=0) only makes sense when acting on functions on S"^ that do not depend on ?/;. Looking at 
the expansion H6.41() of M{x, x'), we see that this is the case when 71 = 0. So can only be allowed to 

act on the terms that have coefHcients ajoo- We gef*^ 



^(„=o)^^(^=o) ^ ^(„=o) J2 amYl,{x)Y'Ux') = amo(Si"=°)YLo(x))YL(^') 



(6.91) 



Im 



From H6.71|l . restoring fi n/k, we see that BxYinni^) = (1 ~ A'^/^^)/Q^iooY^„(a;), and it immediately 
follows that B is indeed the inverse of M{x,x') on this restricted space^^. Before we proceed to the general 
case, let us perform a consistency check: If really is (part of) the Dirichlet-to-Robin operator, it has 

to be first order, and in particular it should have the square root of the laplacian as its principal symbol. 
One (perhaps not very rigorous) way to check this is to examine the limit of (|6.9U|) as Aq becomes large. 
This is easy to do, using the formulas |75j for the asymptotic expansion of the Legendre functions for I large: 



Pi (cosh a) 
So clearly, for I large 



g(i + l/2)Q 

(2A7rsinha)^ 



1 + 



);(coshQ;) 



7r2 e 



-(i + l/2)a 



{21 sinha) 2 



1 + 



2k 



H Pi{k/Mik/ti) 



kl 



l + 0(- 



(6.92) 



(6.93) 



where we used cosh a = A;///, and so sinha = 1/ ^^Jl — /i2/fc2. Thus (noting that Aq ~ P for I large) we 
can conclude that does indeed scale as 



B("=°) = k^o 



l + O 



/A^ 



(6.94) 



''^This is just a formal trick, and could easily be dispensed with. 

"^^This is of course simply the Wronskian of the two linearly independent Legendre functions, given e.g. in 1751 . 

"^^We could have written here ctjmOi but since they are all equal to «;oo we just write the latter. 

'"'The extra factor (1 — fi^ /k^) clearly has to do with the normalisation of M(x,x') which we have left arbitrary. 
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and thus has the right principal symbol to be the Dirichlet-to-Robin operator. 

After this warm-up, we now turn to the n 7^ case. Again we have the sum (I6.84|) of the interior and 
exterior Dirichlet-to-Robin operators, which are given in (|t).30|l and (|(j.35|l . Adding them gives: 



1 — /i-^r^ 



dr In ■ 



Fi{r) 



-^Ao + 1/4 



^3(0 



(6.95) 



4 y l-/i2r2 Fi{r)F:i{r) 
Recall here that we have defined u = (1 + jikr^) / [ijikr^]. To proceed, we use the identity*^ 

2n{l-ii/k) (1 + iikr^\'^^^'^^'' T{2l + 2)T{{1- ^Ji/k)n) 

T{1 + I + n)V{l + I - ^i/k ■ n) 

, we write 
T 



(6.96) 



r(l — /ifcr^) \1 — iikr 
which generalises H6.88|l . Calling T the right-hand side of 



1-^2^2 Fi{r)Fj,{r) 



(6.97) 



Now we finally take the limit r 1/k, scaling A as in H6.90|l to obtain the complete Dirichlet-to-Robin 
operator for QTN: 



B = 



/AT 



r(i + 2VAo + i/4)r((i - ^i/k)^) 



^-i^/kj r(i + ^Ao + 1/4 + VA7)r(-i + ^Ao + 1/4 - fi/k ■ VA7) 



F - - VAo + 1/4 + Va:, - - yAo + 1/4 - f VaT; 1 + (1 - f 



(6.98) 



F Q + \/Ao + 1/4 + VaT, i + ^Ao + 1/4 - 1 + 2VA0 + I/4; 



At this point we really should repeat the check (that we previously performed for ^B^"^'')) of whether the 
principal symbol of B given by (|6.98() really is the square root of the boundary laplacian A b l|A.24|) . Clearly 
in the case Aq ^ Ai H6.98|l reduces to H6.9()|l so we can apply the previous result, but taking other scaling 
limits, like for instance both Aq and Ai large, seems a rather daunting task. Perhaps it is possible to use 
more formal methods to prove that B is first-order, but for the moment we content ourselves with the simple 
consistency check we just discussed, and proceed to the final step. 

The expression for B in (|6.98l) should look familiar. If we act with B on an eigenfunction Yj,j„(x) on the 
Berger sphere, which really means substituting Aq = l{l + 1) and Ai = everywhere, we get precisely*^ 

B,Yl,^{x) = -^Y^„„(x) (6.99) 



with ainn as defined in 
B^M{x,x') = 



So we finally have all the necessary ingredients to show that 



^6.a,„„Y5„„(a;)YL(2:') = E YL(^)yL(^') = ;^^^'^(^'^') 



(6.100) 

— ■ ^ / 

Imn Imn 

which was the main calculational objective of this article. 

Given the discussion in section^] we conclude that M{x, x') (as defined in H2.30|l ') is indeed a conformally 
invariant correlation function on the Berger sphere, corresponding to an operator of conformal dimension 
one. 



Which we can verify by inserting specific values of / and n. 
*^Up to a (1 — p? /k^) term and other possible normalisation factors that do not depend on I and n. 
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7 Conclusions 



In this article we demonstrated that our correlation function M{x,x') (see (|2.30|l ') on the Berger sphere is 
conformally invariant, and is a two-point function of an operator of conformal dimension A = 1. We did 
this by explicitly constructing a conformally covariant pseudodifferential operator, the Dirichlet-to-Robin 
operator B, which we expected to be the inverse of M{x,x') on general grounds, and then verifying that 
this was indeed the case. Note that the methods that led to these two objects were (at least superficially) 
very different from each other: M{x,x') was found using twistor methods |2, namely Penrose's nonlinear 
graviton construction, while B was found by factorising the QTN laplacian, which in the end came down 
to solving a Ricatti-type equation. Thus we believe we have performed an independent verification of the 
conformal invariance of our correlation function. On the other hand, as discussed in section [4. II to obtain 
this result we have made a number of assumptions and generalisations which should eventually be better 
justified. We believe that this will probably require some input from geometric scattering theory. 

We should emphasise here that we are not claiming uniqueness of the correlation function. There could 
very well be other two-point functions on the Berger sphere that correspond to A = 1 scalar operators. 
In fact, the very method we used to prove conformal invariance of M{x,x') also seems to hint that it is 
not unique! Recall [77\ [75] that there are parameter ranges where, apart from QTN, the Berger sphere is 
the conformal infinity of another asymptotically hyperbolic 4-manifold, known as AdS-Taub-Bolt. Now, 
as discussed in (|4.2ll the main feature of the Dirichlet-to-Robin operator is that it is different for non- 
diffeomorphic manifolds. So, although the calculation of the corresponding M{x,x') for AdS-Taub-Bolt 
seems out of reach at the moment, it is likely that it will have a different functional form from the one for 
QTN we have been investigating (although both will have the same conformally flat limit, since M{x,x') is 
unique there). It is intriguing to wonder whether that is all, or whether there are several other conformally 
invariant correlation functions on the Berger sphere. 

Although we focused narrowly on the very particular problem of A = 1 and the Berger sphere, we 
believe that the techniques we employed should be applicable to a wider range of questions that might 
arise in applications of holographic ideas to asymptotically hyperbolic bulk manifolds. A drawback is that 
calculability in our case hinged on the existence and relative simplicity of an eigenfunction expansion, which 
is beyond reach for many manifolds. But perhaps the formal aspects of the Dirichlet-to-Robin operator (and 
related nonlocal boundary operators) viewed as conformally invariant operators will eventually turn out to 
be useful in better understanding conformal field theory, especially on curved manifolds. Of course, as the 
dicussion in section [4.21 tries to make clear, the Dirichlet-to-Robin operator should also play an important 
role in exploring the connections between holography in an AdS/CFT context and the mathematical fields 
of geometric scattering and inverse boundary problems. 

One could conceive looking for the nonlocal boundary conformally invariant operators corresponding to 
correlation functions of higher conformal dimension. Branson and Cover 40 construct higher order, and also 
negative order, generalisations of the Dirichlet-to-Robin operator, which thus have principal parts (A^)^ 
(in (41) this is extended further to arbitrary real order). (We put emphasis on the negative order cases, 
because by the relation 2X ~ d — s H1.4|) , the order of the pseudodifferential operators corresponding to CFT 
operators with the usual conformal dimensions one encounters (i.e. A > d/2) will be negative.) They do 
this by considering conformally invariant boundary value problems for the (bulk) higher order conformally 
invariant powers of the laplacian we briefly discussed in section 13.21 (of which the Paneitz operator H3.7(l is 
an example). There are restrictions to doing this on an arbitrary curved conformal manifold (related to 
the restrictions in defining the powers of the laplacian), but supposing that one can construct conformally 
invariant boundary problems in the bulk, which result in boundary pseudodifferential operators associated 
to a given correlation function, one could ask how this (conformal, nonlocal if the bulk operators are negative 
order) description of the bulk degrees of freedom differs from the usual (non-conformal, local) one based on 
the massive bulk laplacian. 

After these comments, let us return to the Berger sphere. In the end, are we any closer to understanding 
the nature of the field theory that might be dual to gravity on QTN? In an important sense we are, since we 
have demonstrated, in an explicit (but very special) example, that the usual AdS/CFT dictionary between 
bulk fields and boundary operators holds unmodified in our more general case. However, a lot more needs 
to be done to fully understand how holography works for QTN. The best one could wish for, of course, is a 
string/M theory realisation, which should give a better picture of the gauge theory side (of course, given the 
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apparent lack of supersymmetry, such a configuration would be interesting on its own, were it to be stable). 
One could also try to focus on the region close to ^ = k (where the manifold is almost CH^) which we have 
been largely ignoring in this article: Despite the discussions in 13, 79 , holography in the CR case is still a 
bit mysterious and deserves to be better studied. 

It will also be important to obtain correlation functions of higher A, since they might be able to give 
us insight into how the conformal field theory makes the transition from the "stable" region (/z^ < 3/4fc^) 
where the boundary scalar curvature is positive, to the (presumably) unstable region ifi"^ > 3/4fc2) where it 
becomes negative. Recall that one of the main motivations for studying this model is understanding this 
instability (implied by the results of IHOIIS)) and how/whether it is reflected in the geometry of the bulk. As 
an indication of the issues that need to be addressed, we mention a comment in |S2 that, as the boundary 
scalar curvature shifts to negative values, numerically it seems that the bulk laplacian for Quaternionic 
Taub~NUT begins, as expected perhaps, to develop normalisable modes above the Breitenlohner-Freedman 
bound of —9/4. It is plausible that these modes could render the usual AdS/CFT correspondence ambiguous. 
What is more interesting however, is that there seems to be a range in which the boundary scalar curvature 
is negative, but no such modes have yet appeared in the bulk spectrum. Resolving these and other issues 
connected to holography on QTN promises to be a worthwhile endeavor. 
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A More on the Berger sphere 

In this appendix we include a few details on the homogeneously squashed three-sphere, or Berger sphere, 
which we sometimes denote by S^. Its metric is simply 

5,, : [a'ir){al + al) + hirfal] (A.l) 

where the sigmas are the usual left-invariant SU(2) one-forms: 

(Ti = - (cos i/^dfl + sin sin 6'diy9) , 

(72 = -(— sin-i/'df? + cos?/; sin6'diy9) , (A. 2) 

0-3 — ^(dV' + cosOdip) . 

normalised to satisfy dui = —2(72 A (T3, and cyclic permutations thereof. The ranges of 0,Lp and are 
[0,7r), [0, 27r) and [0,47r) respectively. 

The squashing is given by the ratio A(7') = b{r)'^ / a{r)'^ . A(r) = 1 corresponds to the round three-sphere, 
< A < 1 is known as oblate squashing, while A > 1 is called prolate squashing. For definiteness we only 
discuss the prolate case in this article, but our results should hold for the oblate case too, by taking fj, —i^ 
everywhere 

Any constant-r hypersurface in QTN defines a squashed three-sphere, but the squashing depends on the 
radial distance. Also, the overall scale of the metric is important in taking various limits. So we will mostly 
keep a(r) and b{r) arbitrary, and write down special cases at the end. 

The squashed sphere was famously considered by Hitchin |83| in his discussion of the space of harmonic 
spinors (the null space of the Dirac operator) on a manifold: S"^ is a good illustration of the fact that the 
number of harmonic spinors is not a topological invariant of the manifold, but depends on the particular 
metric too. On the physics side, scalar quantum field theory on the homogeneously squashed sphere has been 
considered by several authors |84l 1851 1^ . mainly due to its appearance as (a particular case of) the spatial 
section of the mixmaster cosmological model |HZj- Rather remarkably, the SU(2) x U(l) symmetry means 
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that the problem can be mapped to that of the symmetric quantum-mechanical rotor (see e.g. [SS], chapter 
19, or sections 82 and 103) which (like that of the spherical rotor) is integrable. (As we will see, this 
means that it is relatively trivial to find the eigenfunctions of the laplacian on S^.) There has also been some 
work on S'^ in the context of Kaluza-Klein compactification |^|^ (see also a discussion in ;92'), and the 
asymptotic expansion of the heat kernel for the laplacian has been computed in More recently, Dowker 
and collaborators ,94, have considered the calculation of effective actions for both scalars and fermions 
on S^, aiming at a comparison with the AdS/CFT results of |77[ 178 , 96 . The relation is far from obvious, 
since the regimes where the results are expected to apply are vastly different (see [75] for a discussion). 



A.l Geometry 

Written in terms of the coordinates 0, ip, i/j, the metric (jA.lj) is 



1 

gy =a(r)^( sin^ 61 + A cos^ 6* XcosO | . (A.3) 
A cos 6* A 



The nonzero Christoffell symbols are 



- (A - 1) sin 9 cos 6, = - sin , 

_/ AAcos^ „_AJ_ , . 

^^^-[^ 2y'sin0' ^<'^~ 2sin0' ^^'^^ 

^ _ sin^ 6* + (2 - A) cos^ 6 ^ _ Xcoad 
" 2 sin 6* ' "'^ " 2 shT^ ■ 

The Berger sphere is not an Einstein manifold, but it comes quite close: It satisfies the Einstein- Weyl 
equation 7?,[y] ~ gij, where TZij is computed using the Weyl (non~Levi-Civita) connection. For a recent 
survey of Einstein- Weyl geometry, see On the physics side, it was shown in [SSI that (at least for /i 

close to zero), the spacetime can be self-consistently supported by the stress-tensor of a conformally coupled 
scalar. 

The scalar curvature is given by 

1 4-A(r) , 
a(ry 1 

which for the round sphere (A = 1) of radius 1/fc results in TZh = 3/2k^. Notice that the unit sphere (which 
has TZh — 6) corresponds to choosing k = 2. There are two main cases to consider: If we are interested in 
the squashed sphere at a particular radial distance in QTN, we have 

Then the scalar curvature at a radial distance r is given by 

Tlr = ^;S(r^^(3 - + 4^V^ + ^2,2^4) (A.7) 

This is the expression that is used when factorising the laplacian on QTN in section IHTI If, on the other 
hand, we want to look at the squashed sphere at the boundary at infinity, we have to choose a defining 
function. Choosing (1 — yp'r^) for this purpose, we end up with the metric ljA.l|l with 

a2 = l(l-/.Vfc2), 52 = 1 (A.8) 

then the scalar curvature is 

fcV3-4/iVfc2\ 

^--y( (i-,W j ^^-'^ 
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As mentioned in the main text (and also in [3]), the fact that TZb becomes negative for /i^ > 3/4fc^ should 
have very interesting implications for the gauge theory side. 

Before moving on to discuss the symmetries, we can quickly calculate the Cotton tensor H3.2|l . which, as 
discussed in section^ will tell us whether the space is conformally flat. Of course, it turns out to be nonzero 
apart from the case of the round sphere. 

A. 2 Isometries and conformal isometries 

The squashed three-sphere has 4 Killing vectors, corresponding to its SU(2) x U(l) isomctry group. We use 
the basis ^31 

h = id^ , 

h = i{s\nipde + (cos 6*9^ ~ d^)} ■, 

sin (^-l*^) 

?3 = i{cosipdB ^^(cos6'9y - d^)} , 

sm 

ri = id^ . 

The li satisfy the commutation relations of SU(2), [li,l2\ = ih etc., while ri commutes with all of them. 
The notation comes from the fact that the li correspond to the unbroken (left) part of the bi-invariant 
SU(2)l X SU(2)/j isometry group, and ri is the remaining U(l) subgroup of SU(2)/j. 

For completeness, we include here the Killing vectors for the case of the round sphere, i.e. for the case 
that A = 1 above. Then the four vectors above can be supplemented with two more, 

r2 = j{sinV'9e + '^"^ 1* (cos 9d^ - d^)} , 

sin J ^^-^^^ 

rg = i{cosi)dg —-{cosOd^ - d^)} , 

sm 

restoring the full SU(2)i x SU(2)ii;. Of course, in this case we also have a nontrivial conformal group. The 
conformal Killing vectors can be easily found by restricting the Killing vectors of the AdS4 bulk, and can be 
taken to be 

sin ( ^\ 



ki = 2i sm - cos ( — - — ) de + t g [o^ + o^,) 



2 V 2 



2 V 2 



cos ( ^1 



fc2 = 2z sm - sm ( — - — ) de - i g [o^ + d^,) 



ijj-ip 



.Bin(V 



(A.12) 



^3 = 2i cos - cos ( — ^ — ) de + i — g [d^ ~ d^) 



cos 



V 2 



= 2i cos - sin — - — de - i „ [d^ - d^) 

I \ I ) sm I 

All told, we have ten Killing and conformal Killing vectors, corresponding to the conformal group S0(1, 4), 
which is of course the same as the isometry group of AdS4. We can easily check the conformal Killing 
equation for the fc„s: 

'Cfe„(x)5'ji = Vifcnj + \1 jkn.i = 2ak^{x)gij (A. 13) 

where {x) = V • fc„(x)/3. The coefficient of is of course uniquely determined to be 2, as we can see by 
acting with the inverse metric on (|A.13|I . 

There is another special value of the squashing parameter A: When A — > cxd (/i — > fc) , which corresponds to 
extreme prolate squashing, the three-sphere degenerates to a CR structure, the boundary of the unit sphere 
in C^. This is the boundary at infinity of the Bergman space CH^, which is the coset space SU(2, 1)/U(2). 
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One would expect that, despite the degeneracy of the metric, the Bergman isometry group SU(2, 1) will 
manifest itself as the conformal isometry group of the boundary, as in the AdS case. In fact this was shown 
to be true by the authors of |13|. who explicitly constructed the requisite conformal Killing vectors: 



61 = 2i sin cos ( ) dg + I 2,. e (^v + (2 + cos 9)8^) 



cosf ^1 



62 = 2i sin "- sin ( ] de - i {d^ + {2 + cose)d^) 



sin(V) 



(A.14) 



63 = 2i cos - cos ( — - — ) de + i ) g {d^ - (2 + cos 9)8^) 



cos 



64 = 2i cos - sin — - — dg ~ i ^—^ {d^ - (2 + cos Q)d^) . 

I \ I ) sm I 

The interesting point about these conformal Killing vectors is that they satisfy the conformal Killing equation 
with a different coefficient than the canonical one in IITI : If we again define c!h,X^) — ' ^n)/3, we find 
that 

3 

^b„(x)gij = -^'^b„ix)gij ■ (A. 15) 

As discussed in |13| . there is nothing wrong with this since here there is no inverse metric so it need not 
impose a particular coefficient. 

After this discussion of the two interesting limiting cases of the metric, we will now assume we are at a 
generic value of A. In discussing two-point functions on the squashed sphere, we find it extremely convenient 
to work with variables invariant under the SU(2) x U(l) symmetry group, rather than the Euler angles of 
the two points {9r, fr-, V'r and 9s,ips, ips respectively). We would also, of course, like to work with variables 
that are symmetric under the interchange Xr Xs, since then the symmetry of the two-point function is 
manifest. These two requirements lead us to define the two combinations 

9r 9s 1 9j^ 9 s \ 

Vi = COS — cos — COS -(l/'r + - Vs - ^s) + siu y sin — COS -(■(/'r - -iPs + <Ps) (A. 16) 

and ^ 

V2 = -{1 + cosftr cosf?s + sin^r sin6's cos((y9r — ((Cs)} . (A. 17) 

Any Green's function on can be written as a function of vi and V2 only (six Euler angles minus four 
constraints from SU(2) x U(l) equals two independent variables). One can check that both vi and V2 are 
SU(2) X U(l) invariant, by which we mean that for any Killing vector v G {U, ri}, we have 

However, if we want to impose invariance under the full SU(2)l x SU(2)ij of the round sphere, i.e. allow v 
to be r2 or ra, we find that V2 is not invariant, while vi remains so. This means that on the round sphere 
two-point functions will depend only on vi, as is indeed the case for M'^{x, x') in section |5| 

A. 3 The laplacian 

From the metric (IA.1|I we readily obtain the laplacian 



v-9 



dee + cot 9de + 



- {dipip — 2 cos 9dip^ + d^^p) 





■ 1 1 " 


+ 


52 ~ ^ 



(A.19) 
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Here the subscript B can be taken to stand for "Berger" or "Boundary" . It is useful to define the operators 

1 



-Ao = deo + cot Ode 



sm 



: {diptp - 2 cos 9d^^ + d^^) 



(A.20) 



and 



-Ai = d^^ (A.2I) 

We will also find it convenient to work with a laplacian that is positive definite, and which depends only on 
the squashing (and not on the overall scale). So we define an operator A(,.) as 



(A.22) 



Given the definitions of Aq and Ai above, at the inner boundary located at a particular radial distance r, 
we have (recall (|A.6|l 'l 



A(r) = Ao + 



2 22 - - k"^', 
1 ~ fi^k^r"^ 



-Ai 



At the boundary at infinity, given by r ^ I/fc, we have, from (|A.8p or directly from l|A.23|l 

,2 



lim A 



(r) 



(A.23) 



(A.24) 



A. 4 Separation of variables 

To look for eigenfunctions of the laplacian (|A.I9|I , we separate variables in a standard way: 

YL(0,^,V') = e-''"'^e-™^UW 
with "rfmni^) to be determined. The variables f and '0 separate immediately and we are left with 



(A.25) 



do 



cot Ode 



sm 



-{rn? — 2mncos9 + n^) 



1 

62 



Yl„ix), (A.26) 



where x stands for {6,(p,ip). Renaming cos( 



we obtain 



1 



(1 - z')d. 



2zd, 



m — 2mnz + n 



1 



1 

62 



(A.27) 



The first bracket is a well-known differential equation, arising from representation theory on SU(2), and its 
eigenfunctions are the Wigner functions ( 88 , see also [^|73]), sometimes denoted d'^„(z). We follow [7^1 
in denoting them P^^^{z), a notation that indicates that they are simply a generalisation of the Legendre 
functions P;(z) and associated Legendre functions P'„(z). The eigenvalues of the Wigner functions are simply 
+ 1), depending only on the principal quantum number I. In other words, the first bracket of ljA.27p 
acts on P^^niz) as: 



(1 - z^)d,, - 2zd, - 



m — 2mnz + n 

r^"^2 



PLn{z) = -l{l + l)PLn{z) • 



(A.28) 



We can call this equation the "SU(2) differential equation" since it always arises when doing representation 
theory on SU(2), as for instance when describing the spherical top. As discussed in the beginning of this 
appendix, it is a nice (and well-known) result that the Wigner functions are also eigenfunctions of the 
more general equation (|A.27() corresponding to the homogeneously squashed sphere, or equivalently to the 
symmetric top, with very simple eigenvalues^''. So substituting rf^-^^^{z) — Plnni^) i'^ l|A.25|) . we have that 
the Yj„„(6',(^,'0) satisfy 



I 

62 



n 



YL(^^,^,V') 



(A.29) 



^"^This is no longer true when we consider the asymmetric top, or equivalently the full mixmaster solution a^af + b^(T| +0^17 
There the analysis becomes much more involved and can only be approached numerically |84| . 



3„2 
3- 
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As mentioned earlier, we find it convenient to use the positive laplacian A(r) defined in IIA.221) rather that 
V^. The operators Aq and Ai we defined above clearly give 

AoYL„{9, V, V-) = l{l + 1)yL(0, V, ^) (A.30) 

and 

AiYj„„(0, ^, V-) = ii?Yl„{B, ip, V) (A.31) 

so that 

A(,)Y,L(e,¥',V') - + 1) - ^''''fi I ''^ } V. ^) ■ (A.32) 

Since A(r) is finite as r 1/fc, if we want to look at the Berger sphere at the boundary at infinity (as in the 
eigenfunction expansions of sections [S] and (HI , we can just take r = 1/fc in A^. We thus obtain^^ : 

AbYL„(0,¥>,V) = + ^ •n2|Y5„„(0,^,V) . (A.33) 

If the range of / is taken to be ^ = 0, 1, . . . oo, and that of m and n as, —1,-1 + 1, ... ,1—1,1, the Yj^j„(0, Lp, ijj) 
form a complete orthonormal basis on S"^. They satisfy the following orthogonality relation 

7T /•27r /•47r ^ IQn^ 

Y'i^i„,(6',<^,V')Y™^„j(6',93,V')d-0d.^sine'd6' = ^^--pY5;ii2(5™i™2(5„i„2 , (A.34) 



IQ JO JO 

and also completeness 



E T6^Y,L(0, V, m\.n{o', v': - 5{e - e')5{^ - ^')d{^ - r) . (A.35) 

Imn 

For this reason in the main text we actually use a normalised version: 

yL„(0,^, V) = y2me-""*e-™'^Pi„(cos0) . (A.36) 



B Pseudodifferential Operators on Manifolds 

Here we will give just a few basic definitions from the theory of pseudodifferential operators. The reader 
interested in better understanding pseudodifferential operators should consult one of the books on the subject 
(e.g. ,98, 99, 46 ). See also 100^ for a discussion oriented towards applications to index theorems. 

The definition of a pseudodifferential operator on IR'' usually starts from the observation that one can 
think of a differential operator P acting on a function f{x) as defined through its Fourier transform: 

Pf^^) = Trvil ^"M^^Ofim- (B.i) 

JW^ 

Here / is the Fourier transform of /, and the function p{x, ^) is called the symbol of P. Now for differential 
operators, which for order m can be written as P = 'Y^™^Qai{x){—idxy , the symbol is a polynomial in ^: 
p{x, = YlT=o ■ It is natural to generalise by considering symbols that are not polynomials, and this 

leads us to pseudodifferential operators. In that case one has to be more precise about the definition of the 
symbol of the operator, and thus its order. We will not need the full definition, but only a more restricted 
one called a classical symbol: We say that a (classical) pseudodifferential operator belongs to §™ (the set 
of pseudodifferential operators of order m) if its symbol p{x, ^) is smooth in x and ^, and is an asymptotic 
expansion in homogeneous terms, the degree of each term decreasing by an integer. That is, 

oo 

p{x,i)=Y,Pm^,{x,0 (B.2) 

1=0 

®^Note that this is still positive since /j. < k and for unitary representations n < I. 
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where Pm-i is a homogeneous function of degree m — i. This definition is sufficient for our purposes, since, 
as discussed in the Dirichlet-to-Neumann operator is a classical pseudodifferential operator. 

Another useful definition is that of the principal symbol of a pseudodifferential operator, which is the 
highest-order part pm of its symbol. The principal symbol contains a large amount of information about 
the operator, for instance it is enough to decide whether the operator is elliptic or not. 

In the theory of pseudodifferential operators, one often encounters infinitely smoothing operators, the 
operators whose symbols belong to S~°° (defined as the intersection of all the S"*). These operators will 
completely smoothen out the support of any distribution. They can be arbitrarily added to any asymptotic 
series (like (jB.2|l ). so one is often forced to work modulo smoothing operators. 

Probably the major conceptual difference between differential and pseudodifferential operators is that 
the latter are nonlocal. However they are a very restricted set of nonlocal operators, being what is known as 
pseudolocal. This means that, although they are defined through the Fourier transform, which smears out 
the support of the function they act on, they don't smear out the singular support. In other words, if a 
function is smooth in some domain, then P acting on it will also be smooth in that domain. 

It is a well-known result, due to Seeley, that any complex power of an elliptic operator is pseudodif- 
ferential. This tells us that the operator (|5.9|l is well defined. As for the expression we found for the full 
Dirichlet-to-Robin operator on QTN, we will need the following result (see e.g. p. 297): A function 

p(Ai, A2) of two self-adjoint operators is pseudodifferential if Ai and A2 (in our case a/Aq + 1/4 and 
commute, Af + A\ is elliptic, and the function p{X) is a symbol in the appropriate symbol class. The first 
two conditions are trivial to show, however we have actually shown the third condition (cf. (|6.94|l ') only for 
a special case. 

So far we have been discussing pseudodifferential operators on H''. However they can be readily defined 
on curved manifolds, by first defining them on a coordinate patch (transferring back to K''), and then 
requiring that the definition holds for all the charts in an atlas. We will not require much knowledge of 
the technical issues, however it is worth noting that only the principal symbol is diffeomorphism invariant. 
Clearly this makes it rather tedious to work with pseudodifferential operators in a covariant manner, and 
there has been a lot of effort to try to overcome this difficulty (see e.g. jlUl) for a discussion). 

Fortunately in our case things are rather simple: As dicussed in appendix^ the Berger sphere allows an 
explicit expansion into known eigenfunctions, and the eigenvalues are simple enough that we can directly act 
on these eigenfunctions with any pseudodifferential operator by naively applying the spectral theorem. That 
means that, assuming we know the eigenfunctions of some (pseudo) differential operator A, we write 

AY'' = \aY-' ^ piA)^ ^ p{\a)Y'' . (B.3) 

Apart from the very important application of pseudodifferential operators in the Atiyah-Singer index theorem 
(see e.g. [lOOp . they have long been useful in quantum field theory in curved spacetimes (see e.g. |102| '): For P 
an elliptic, positive order operator, the heat kernel e~*^ and the zeta function Tr{P~^) are pseudodifferential. 
As for the square root of the laplacian, it has also appeared in quantum field theory in various contexts, mainly 
related to the construction of non-local actions (some sample references in this field are |lU3lll(J4lllU5| '). | 1(J4| 
actually discusses a factorisation of the (fiat) laplacian similar to that in section ITTl In its lorentzian, fiat 
space edition, the square root of the wave operator in 2+1 dimensions is known to satisfy Huygens' principle 
(unlike the wave operator itself) |1U6I I1U7| and thus it is natural for it to appear in 2+1-dimensional 
bosonization (e.g. |ll)8j ') (and thus also invalidating to some extent our comment in a footnote in the 
introduction that the "square root of the laplacian is nonlocal so has nothing to do with the Dirac operator"). 
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